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■ CDMA Patent 



The process at the transmitter of CDMA- is roughly 

BitStream -> Frame Builder — >GoldCode Encoding -> Modulator [Using LookUp Table] 

D/A -> Fihcrf Switched Capacitor, Clock Frequency] -> Fi!tcr[ ResistorSjCapasitors] — > Transmitter 

The corresponding process for CDMA receiver is 

Complex Number Seq -> CDMA Process -> Channel Estimation -> Demodulation [Using Viicrbi] -> BitStream 
These process arc in fact similar to GSM so wc could build a dual mode phone. We describe the GSM processes 
Tlie transmitter processes are listed below 

BitStream — > Frame Builder — > Modulator [Using LookUp Table] 

D/A -> Filler[ Switched Capacitor.Clock Frequency] — > Filter[ Rcsistors,Capasitors] — > Transmitter 



At the receiver wc have the reverse process 

Complex Number Seq -> Multiplying incomming stream by i'^n -> BitSinchronisation -> Channel Estimation -> 
Demodulation [Using Viterbi] -> BitStream 

Note 

1) TIie mode dependent process are highlighted. 

2) The same filter can be used in both the cases if the clock frequencies are selected with care. 

3) Thc liming issues, and frequency dependent items like antennas have been ignored in this paper. 

4) ln this paper we use many different pulses depending on the bandwidth, chip rate, required and we trade off the non 
constant amplitude. This is the key inovation (to get higher bit rates). 

The selection of bitrate is just a matter of engineering design involving BER, Power Amplifier efficiency and spectrum 
utilisation. Wc also point out that although in this paper wc have designed a modulator from first principles, wc would 
in pralitc use the current Nokia Mobile Phones lookup table archilccurc 14T or 6TJ obtaining the correct syclostalioary 
spectrum in each mode by the appropialc filler [Analog]. The Switched Capacitor will in praiicc do most of the shaping 
common to both. This again is subject to design 

5) Thc key inovalivc step is the nature of the function used to iransform the gold code into the sequence used in the 
correlator by the function below. It is this function that we Vkish to protect in the patent. The transformation need only 
occur once.and the results stored, as the mobile can pre calculate the sequence di for i = 0.1. .. N-1 given die code 

{ cq , C] , Cfi _i } when the code is asignmcnl. 



yi = (-1)' fori = 0. 1. 2. .... N - 1 

Given code { Cq, Ci cn-i } where N is the number of elements in the sequence . 

ai .= 1 if Ci = I; fori = 0, ...N - 1 
aj = -1 if Ci = O; fori = O, ... N - 1 
bo = Bo : 

bi = bi_i + a. fori = 1,2 N - 1 : 

d; = yj i"' for i =i 0. 1 , 2 N - 1 and 

f = V(-l) 
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In, e function is behaviour is ispecified in Ihe variable x4 in the CDMACoarseSynchroniser. 

We also stress that while the biLstream is differentially encoded and ihc subsutuiion used is { O -» 1 ,1 -1 } the 
substitution elsewhere used for the gold codes is {0 ^ — 1 , 1 -» 1 } 

6) The demodulation in the CDMA mode docs not need a Vilcrbi alfiorithm as is tequircd by the GSM mode. Thus most 
of the work is done by the gold code, during correlation. 

7) There is nothing special about the Gold Code used in the example. Any gold like code will do as long as it possescs 
the usual properties of a gold code. 



nnr_nnr I nfc_il 7flJ n7JJJ0 QJ71 Ofl PC:qi BR-Cfl-Cfl 



CDMA CSMIn troductionPutent. nh 



/The funclion is behaviour is specified in the variable x4 in the CDMACoarseSynchroniser. 
We also stress that while the bitsiream is differentially encoded and the subsUlution used is { O -»1 .1 -» -1 } the 
substitution elsewhere used for the gold codes is {0 — 1, 1 -» 1 ) 

6) The demodulation in the CDMA mode does not need a Vitcrbi algorithm as is tequircd by the GSM mode. Thus, mi 
of the work is done by the gold code, during correlation. 

7) Thcrc is nothing special about the Cold Code used in the example. Any gold like code will do as long as it posseses 
the usual properties of a gold code. 
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■ Dual Mode CDMA and GSM Operation 



In this work we demonstrate how to build a dual mode phone using as much as possible the same haidware for both the 
systems. 

Given the following process for GSM Transmitter 



tream — > Frame Builder -> Modulator [Using LookTJp Table] 

D/A -> Filtef[ Switched Capacitor.Clock Frequency] -> Filter[ Resistors,Capasitors] -> Transmitter 

At the receiver we have the reverse process 

Complex Number Seq -> Channel Estimation -> Demodulation [Using Viterbi] -> BitStream 
The corresponding process at the CDMA transmitter is roughly 

BitStream -> Frame Builder ->GoldCode Encoding -> Modulator [Using X^okUp Tabic] 

D/A -> Fillcr[ Switched C^aci tor, Clock Frequency] -> Filter[ Rcsistors.Capasitors] -> Transmitter 

The corresponding process for CDMA receiver is 

Complex Number Scq -> CDMA Process -> Channel Estimation Demodulation [Using Binary Decision] -> 



Note 

1) Thc mode dependent process arc highlighted. 

2) The same filter can be used in both the cases if the clock frequencies are selected with care. 

3) The timing issues, and fr^equency dependent items like antennas have been ignored in this paper. 

4) In this paper we use many different pulses depending on the bandwiddi, chip rate, required and we trade off the non 
constant amplitude. This is the key inovation (to get higher bit rates). 

The selection of bitrate is just a matter of engineering design involving HER. Power Amplifier efficiency and spcclnim 
utilisation. Wc also point out that although in this paper we have designed a modulator fi-om first principles, we would 
in pratice use the current Nokia Mobile Phones lookup table architecure [ based on a pulse width of 4T ] or an 
architecure based on 6T pulse width obtaining the correct syclostatioary spectrum in each mode by the appropiate filter 
[Analog]. The Switched Capacitor will in pratice do most of the spectral shaping common to both modes. This again is 
standard to design 

5) The key inovative step is the nature of the function used to u-ansfonn the gold code into the sequence used in the 
correlator by the function below. It is this function that we wish to protect in the patent. The transformation need only 
occur once.and the results stored, as the mobile can pre calculate the sequence di for i = 0,1. N-1 given the code 
(co.ci, cn_i ) when the code is asignxncnt. 



yi = (-1)' fori = 0, 1, 2, .... N - 1 

Given code ( co, C| , Cn-i 1 where N is the number of elements in ihc sequence. 

ai = lif Ci = 1: fori = 0, ...N-1 
aj = -1 if Ci = 0; fori = 0, ... N - 1 
bo = ao ; 

bi = bi_, + Eiforl = 1, 2 N - 1 ; 

di = yi fori = 0, 1, 2 N- 1 and 

i = V(-i) 
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-J^l]:= nBBds [■■RixiaP*uxic;t:i.oxxBl]l££Em:odBd'' ■■ ] 

ihlSJ Kain«B ["RinBrFuiictl.oiieDl££Er»cod««a^*"] 

Outf5i= {AllGoldse<3uences, Autocorrel^tionSequence, 

AutamorpliismSigiua, CrosscorrelationSequance, CyclicMulciplativeGroup, 
DropLeadingZeros, Gold^equence, GraefffeMethod. InicialConditions , MinimvunPoly , 
ModuloMuatplication, PolynomialKultiplication, PoasifaleDi visors. Ring-Division, 
RingPower, SequenceGenerator, Specif iedGoldSei^uences, TraceRepr-esentacion, 
TupleRepresentacion, TwoAdicExpansion, UnitsRing, ZeroPad, Z^s-oSequencee , T. a) 

xnl3] := Needs ["I.a\ii:ontPuaiccloiifl'"] 



ated when Needs was evaliiaced. 



Tnf4J ;= J»Be<ls["r.awreatfcioc»tionTeBt'"] 

Needs: ;noconc. : Context Lavureni^JNotationTeet • was i 

Information on ihe functions used can be obtained using help. 
Zn{5j:= Names [■r.auz-entruzxctii.oi^^*"} 

OutCSj= {AKN, AlphaKI. ANKlnitialStateSetUp, BT, FiltPulse, h, hFiltered, Initialstate, J, 
Laurentc, Ls-wz-entLK. Laurents, M. ModulatingPulse. Modulacionlndex. Modulator, 
NumbexOf Curves , Pha^eAngle. PhaseAngleFas t. Receiver, ReceiverProper , S, 
Samplinglnterval, StartingQuadrant, SyncSample, T, C, 4, i(r} 



In (6]:- UodulBtloxxZndex := — • 
Xn't9}:= « ModtilBtiox'Oat:a..in; 

I^fIi;.-= PlotCOptPulseCI,] [0] [t:] , {t, D, 
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The unit of time is T = 1 for OptPulscWe scale the Pulses to T = -axAoF ""^^ of lime. 

Xn[12J := OptPua.BeSc«l«fl[8] [D] [t_] ! = OptPixlso [t] [0] [t / T] 

OptPuisoBcfti^are] [1] [e_] r= optPui8o[i.] fi] [t /t] 

T/)fI4J.-=, S-lotlOptPulseScalwdCl.] [O] [t] . {t, 0. B T) ] 




5^3^0"®- ^°^^°°°°^2°*^0®20000a 500003 
auc[14]= - Grapliics - 



IntlS} := 
Outn5J = 



CI. 1, -1, -1, 
1. -1, 1, -1. 
1. -1, 1, -1, : 
1. 1, 1, -1, 1 



1, 1, 1, -1. 1. -1, 1, -1. 1, -1, -1. -a, _i, 1, 1, 1, _i ^1 i -, 

1, 1, -1, -1, 1, 1. 1, 1, -1, 1, -1. 1, _i, 1, _i, 1, 1, 2.. 1, -1. ' 

. -1, 1, 1, -1. -1. 1. -1, 1, -1, 1. 1, 1. 1, _i, _i, 1^ _i 

-a, -1, 1, X, 1, -1. -1, 1, 1, 1, 1, 1. 1, 1. 1, 1, -1, _i, 1, 1 _ij 



laflGJ:^ Plot[FiltPulBe[I.] [O] [t], {t, 0, 8 T) ] 




g^j^Q-©. OOOHDlOOOOaSOOOQD 2000012 500003 
ti[i6j= - Graphics - 



■ The Gold Sequence Set 

We generate the sequences using the method spccilicd by SenJw Bortas and P Vijay Kumar in Ref [1]. The numbering 
of the sequences is the one used in the paper. We generate a small subset of the sequences. There arc 2^°+ 1 sequences 
with the quaternary polynomial used. Given any binary primitive polynomial, we can generate the corresponding 
quaternary polynomial. 
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J Sp^cl£iedGolas«q:uexicao[(l, 3, 2, 1. 0, 3, 0, 0, 1)] 2, 3, Laat -1. t^ot}]; 

Thc last sequence of seqiist has nice autocorelation properties. 
The last sequence is in fact a m-sequcnce of lenglh 1023 bits. 
Jiit23] := GoldBetjliet: // uaat 



{1. 


0, 


0, 


0, 


0 . 


0 , 


0. 


0. 


1, 


1, 


0. 


1, 


1, 


0. 


1. 


0. 


1 , 


1 , 


1, 


0 


1, 


0, 


1, 


1. 


0. 


0, 


1, 


0, 


0, 


1, 


1, 


1, 


1, 


1. 


0. 


1. 


1, 


0, 


X. 


1, 


0. 


1, 


0, 


1, 


a. 


0, 


1, 


0, 


1, 


0, 


0. 


1, 


a. 


0, 


0, 


0, 


0, 


1, 


1, 


0, 


1, 


1. 


0, 


D, 


1. 


1, 


1, 


1, 


1, 


0, 


1, 


1, 


1, 


0. 


0, 


0. 


1, 


1. 


0, 


0, 


1, 


a. 


0. 


1, 


1, 


1, 


1, 


0, 


1, 


0. 



1, 0, 1, 1, 0, 
0, 0, 0. 1, 0, 

0, 0. 1. 1. a, 

1, 0. 1. 0, 0. 
0, 1, 0. 1, 1. 
0, 1.0. 1, 0. 
0. 0. 0. 0, 1, 

0, 1, 0, 0, 0, 

X. 0. 1, 0, 1, 



0, 0, 0, 0, 0, 1, 0. 0, 1, 1, \ 

1, 1, 1, 1, 0, 0. 1, 0, 0 0 
0, 1. 0, 1, 1, □. 1, 1, 0, 1, ( 

0. 1, 1, 0, 1, 1, 0, 0, 1. 1. : 

" " ■ ' ' 1. 1, 



0, 0, 0. 

1, 0, 0, 0, I 

0, 0, X. 0, : 
0, 0. 0. 1. : 



1- 0, 



OMt[2i]= {255. -1, -1. 



-1. -1, -1. 
-1. -1, -1, 
-1. -1. -a, ■ 
-1, -1. -1, 



-1, -1. -1, -X. -1, . 



-1, -1. -1, -1, -1, -1, _i, -x. -1. 

-1. -r. -1. -1. -1, -1, -1, -1, -1. -1, 

-1. -1. -1. -1. -1, -1. -1, -1, . 

-1, -1. -1, -1. -1. -1, - 



-I, -1. -1, -1, -1, -1, . 



. -1, -1, -1, -1, 
-1, -1, -1, -1, -1, -1. -1, 



-1, -1, -1. -1, -1 



-1, -1. • 
-1, -1. 

'1- -X. 
-X. -1. - 



The ihird sequence in ihe list has the following autocorrelation 
Xnl25J := Aut:occ>x'3r«lat:dLorLSec2U0zice[GoldseqllsC // # [ [3] ] &] 

Out[25}= {255. -1. -1. -1, -1, -1, -1, -17.. -1, -1, -1. -17, 
-1. -1, 15. -1. -1. 15, 15, -17, -17, -1, -17, 15, - 
15, -1. -1, 15, -17, -1, 15, -1. -1, -17, -17. -17, -1, - 
-17, 15, -1. -1, -1, -1, 15. 15, -1. 15. -1, -1, -17, -i: 
-17, -1, -1. -17, -1, -1, 15. -1. 15, -17. -1, -1, -17, - 
15, -1, 15. 15, -1. 15, -17, 15, 15, -1, 15, -1, 15, -17, 
-1. -17, -1, -1. -1, -1. -1, -1, -17, -1, 15, -17, -17, 1., 

15, -1, 15, IS, -17. 15, -1. 15, 15, -1. 15, -1. 15, -17, -1, 15, -1 -i - 
-1. -17. 15, -1. 15, -1, -1, -17, -1, -1. -17, -1, -1, 15, 15 -1 -i 15 
-17. -17. -1. -1, 15, -1, 15. 15. -1, -1; -1. -1. 15, -17 if ^1 -X if' 
15, -17, -1. -17. -17, -17, -1. -1, 15, -1, .17. 15^ .^^ 15 £5 

-1. -1, 15, -17. -1, -17, -17. 15. 15. -1, -1, 15. -1. 
-17. -17, -17, -1, -17. -1. -1, -1, -17, -1. -1. -1, 



-17. -17, -17, -1, -17 
-1. -17, 15, -1, 15, -1 
-17, 15, 15, 15, -1, - 
15, -1. -1, 15, 15, 
, -1, 15. -1, -17. 1 
L, -1, -1. IS, -IT. 

17, 1 



,15, 
-1. -1, 
i, -1, 
•17. 15, 



15. -1. 15, 



Wc generate the output of the modulator 



In[93}:= ModlOucpuc = MoctuXatojrtl.) [ (GoldB.qliflt // taat) /. {0 -> -1), NUnibBrOf Cur 
Modulat:lxigPulBe Fl-lbPulae, SampljLXL0Znt.ez-iral -* T / A] ; 

We check the output 
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Now we iry to test the modulated sequence using the 
The number of samples per chip is equal to 



3 

26000DOO 




32 

Sei^pli^^Xn'eearira.l -» T/«, KQaultitizagPulBe -r FlltPulae] //Dx-op[#, - 



( ((3alclso«3li»t // i 


asl 


: // Dirop[#, 


/ 


. { 0 - 


> 


-1}) 
















{0, 0. 0, 0, 0, 0, 0, 


0. 


0, 0, 0, 0. 


0, 0. D, 


0, 0. 


0 


0. 0, 


0. 


0, 


D, 


0, 


0, 


0, 


0, 


0, 0, 0, 0, 0, 0. 0, 


0. 


0, 0, 0, 0, 


0, 0, 0. 


0. 0. 


0. 


0, 0. 


0. 


0, 


0, 


0. 


0, 


0, 


0, 


0, 0, 0, 0, 0.0. 0, 


0. 


D, 0, 0, 0. 


0, 0. 0, 


0, 0, 


0. 


0. 0, 


0, 


0. 


0, 


0, 


0, 


0. 


0, 


0, 0, 0, 0, 0. 0, 0, 


0, 


0, 0, 0, 0, 


0, 0, 0, 


0. 0. 


0 


0. 0, 


0, 


Or 


0. 


0. 


0, 


0. 


0. 


0, 0, 0, 0, 0, 0. o. 


0, 


0, 0, 0, 0, 


0, 0, 0. 


0, 0. 


0 


0, 0, 


0. 


0, 


0, 


0 . 


0, 


0, 


0. 


0, 0, 0, 0, 0, 0, 0. 


0. 


0, 0, 0, 0, 


0, 0. 0. 


0, 0, 


0 


0. 0, 


0. 


0, 


0, 


0, 


0, 


0, 


0, 


0, 0, 0, 0. 0. 0, 0, 


0. 


0, 0, 0, 0, 


0, 0, 0. 


0, 0, 


0. 


0. 0, 


0, 


0. 


0, 


0, 


0, 


0. 


0, 


0, 0. 0, 0, 0. 0. 0, 


0. 


0, 0, 0, 0, 


0. 0, 0, 


0, 0, 


0. 


0, 0, 


0. 


0, 


0. 


0, 


0. 


0, 


0, 


0, 0, 0, 0, 0, 0. 0. 


0, 


0, 0, 0. 0, 


0, 0, 0, 


0, 0, 


0. 


0, 0, 


0. 


0, 


0, 


0. 


0, 


0. 


0- 



We have successfully demodulated ihc bilstreani 

inC95j!=^ Mesaout:priitOpt: = Hoctula^or [Jj] [ (ooldBeqlxsc // z.&at) /. (Q -> -1}, 
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ovtzlsei^ - Graphics - 



{H«ceiinar[I,] [Mo<aoucputc>pt, StarClxxffQuadrazib -> 0, 
( (Ooldseallet // Last // »i-oi>[#, -4l6) /. < O -> -1}) 

0, 0, O, 0, 0, 0. 0. 0, 0, 0, 0. 0. 0, 0. 0, 0, 0, 0, 0, 0, 0 0 0 0 



0, 0 

0, 0, 0, 0, 0. 0. 0, 



0. 0. 0. 0, 0, 0, 0. 0, 0. 0, 0, 0, 0, 0, 0, 0, o! o' 
"""""""""""--O, 0, 0, 0, o' 



0, 0.>0, 0, 0, 0, D, 



0, 0, 0, 0. 0. 0, 0. 0. 0, O, 0, O. 0. 0, 0. 0, 0, 0, o' o' o' o' o' o' 

0. 0, 0, O. 0, 0, 0. 0, 0, O, 0, 0. 0, 0, O. 0. 0, 0. 0, 0 o' d' o' o' o' o' 

0, 0, 0, 0, 0. 0, 0. 0. 0, 0, 0. 0. 0. 0, 0, D, 0, 0, 0, 0. 0 o' o' c' o' o' 

0, 0, 0, 0, 0, 0, 0, 0. 0, 0. 0, 0. 0, C3, 0. 0, 0. 0, 0, 0 d1 o' o' o' 0 0 

0, 0, 0, 0. 0, 0, 0. 0. 0. 0, 0, 0. 0, 0, 0. 0, 0, 0, 0, 0. 0, 0 oi o' o' o' 0 



■ CDMA OPERATJON 



As an example given a symbol slream e.g. _i _i, j consisting of -1. and I. the following function 

demonstrates the encoding process. 

xn ( 49] : = csHASncQde [BiPolaia:BitSBQ_r OoiaBe«_] : = 

3tX = GQleSOCT/. {0 -+ -1); 

»a«p[xl#&, BiPolairBltSoii] // yXatten] 

In[34] i= cnMAEncoawdSoti = CDWAEncode [ { - 1 , 1, l, -i}, Ooldseqliat // i.b«c] r 



■ CDMA decoding of single Symbol 

The modular output associated witii ( 1 ) 

Ta[50}:= MoaoutjmtPXuBOne = MoflulatQX- [L] [ CDMnEncode [ {1} , (Goldaatzllst: // X.aat) /. {0 -> -X} 
NUmberOfCuz-ves -» 2. Moaulat^lngtmlBe -» OpCPuXaeScaled, Saa^lijigXzit^ex-v^l -•' T/4]; 

This is a primitive decoder built to study the autocorrelation. This will help in decoding 



SBB-qor 66E-S .Old 
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-^151} Pr-iin±tiiroCI>MlkRecoi-»r«ir-[MoaOii.tiiro.t_ , Oo3.aseiz_, SBiDple_, OvexSainplJjaffl t = 

3cl = Partition [Moaoucput, OverSanvldLntr] // Tranapooo //# r [Sample] ] &» 
3C2Bt:Bte = Goias»q[ / . {0 -> ■* 

x5Bt«t© = TaJbl^r C-l> "Modti, 2], (i, o, Longtlirxastato] -1)]; 
3c3llpdato i= JSodule [ { } , 

x2Bt:at:e = HatateRi^liC C9c2 State ]; x4 = rol<ll.iat [Plus, 0, 3e2scate) //n©8t//l« 
xSBtatQ = RotateRight [xBStatB] ; xl (xSState x4) // Ai.i>ly[pluo, #]&],- 
Tal>le [xBUpdaCa, (1, 1, l.eiigtii[aol«aso«j] } ] ] 



Jn[52J:= Tern = Pz-dLiiiitltreCXaiAn*c«ive^CHotiOuts>utFlueO&e, (GoldseqliBt: // l.ast) , 
Xnf53j:= XilatPlot [Tom // xm, FlocJoiaea -> True. PlotRange -> All] 



1, 4].- 




OutfS47= - Graphics - 

Xnr5S]:= IjistPlot [fDom // Im // Taka [#, 20]6, PlotJoirted -> •. 



. PXetRaaae -> All] 



150 




12 3 








7B 




A 35 




2 . 5 


7.5 10 iTrs 15 i^-5^^0 



Out [55] = - Graphics - 
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= Take [Tom, 10] 

{1.75501 I- 0.329995 I, 0 . 03«3 42 1 - 1 . 583 66 I. 1.39349 + 29.6902 1, 153.017 - 0 95706B 3 
-0.71B377 + 1S5.563 I. -43.6424 - 2.04777 1, 2 . 2103 1 + 0 . 2 5663 6 1 , 5.392+2 16642 1 
-3-13017 + 6.54167 1, 3.61193 - 2.5452 1} , . J:>-<: + . AOO^^ X . 

I-i^tPlot [Tom // Re // Taifco[#, 20]&, PlotJoi-ned -> Tinio, PloeRan0s -> Aii] 




10 12.5 ~ 15 TT-rS' 20 
ovtlSi]" - Graphics - 

Xnf5fl;-= tlottlot [Tom // Xm // T'»ko[#, 20]K, FloC Joitted -> True. PlotAango -> All] 

150jl 




7.5 10 12.5 15 17.5 20 



We try a less favourable sequence 

Jnfffo; .-= MoclOutputPiuBOnea = Modulator [I.] [ CBHABncode [ (1} , t3oldB««3li.st // « [ [3] ]&] , 

NvanbexOfCuz-ves -> 2, Hodulatln^Pulse -» Opt:PulaeScaled« SaapllxiBZzifcBzval -* t/4]; 
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}:= Tom* = Prim±ti-ireCDItuuieceiTrB3r[Mo<30ut;putPluaOiiB3, (Goldawqlist // # [ [3] J t) . i, d] ; 
rnZ-es;— Take [TomA, 10] 

cfff^/= {14.7675 - 4.24542 1, -0.301874 + 8.61806 1, -5.92939 + 29.6038 1. -154.61 - 10 d754 I 
6.6413B + 166.152 I, 45.-2887 - 7 . 16233 I . - 12 . B09 + 2 . 47796 1 , 3 4S94B - 19 7 653 1 
14.2097+6.33624 1,0.565872-6.10244 1} ,J-4^S4B 19.76531. 

(70J.= LiBtPlot[Toin4 // He // TB.ke[#, 20]t, PlotJoinaii -> Tmo. PlotRanfft. -> All] 



-ito|- 

Out(70j= - Graphics - 

Tji/7i;.= LlatPlot [Tom* // Im // Take [#, 20]{i, tlotaoiaed -> Xania, Pic 

15 Oj 




'X^/^itl 3:2:^515 17Sj^ 2 0 



■ CDMA Decoding of Several Symbols with Training sequence Receiver 

First we generate ihc modulator output. For simplicity we will use r very short training sequence . Let the training 
sequence be(l -1 ,- J . 1 J . JLet llie guard sequences bc { 1, 1, 1 ) . JLci the data symbols be generated by a random 

Znm}:= clat:al = {1. 1, X. 1, 1, 1, 0, 1. 0, 1, o, 1, 1, Ir O, 1. 1, O, 3-, 0} 

Inr737:= data2 =. {0, 0. 1, 0. 1, 1, 0, 1, O, T, 1, O, X, 1, D, i, O, 0, X, X)/ 

Iti[74]!= gruar-d = {1. 1, 1) 

Oatf7tfJ= {1, 1, 1} 

I/3f75;.-= trai-niaff = {0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 0, O, 0, 0, 1, O. 0, 0, 1, 0, 0, 1, 0, 1, 1, 1}; 

The GSM Training sequence is ( 0,0. 1,0.0, 1,0, 1,1, 1,0.0.0,0, 1.0.0,0,1, 0.0, 1,0.,1, 1. 1 ). In fact any short m-sequence can be 
used to characterise the output. 



no d 
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-~^ly at ihis point that we difJcrentially encode using the gsin scheme 



In [76]-.= fzri 



•ToinEguaLard, datiBl, rz-alixlng', aa.t:B2, guard] 



IB error: new symbol name "frame" is similar 

OuC[76J= {1, 1, 1, 1, 1, 1, 1. 1, 1. 0, 1, 0. 1, 0, 1, 1. 1. 0, 1, 1 

0, 1, 0, 1, 1, 1, 0. 0, 0, 0,1, 0, 0, 0. 1, 0, 0, 1, 0 l' 

1, 0, 1, 0. 1, 1, 0. 1, 1, 0, 1, 0, 0, 1. 1. 1, 1, 1} ' ' 

Jn[ei} ;= OSKD±££E^c:cicleaFz'aiiie [ farM»e_] : = Hodule [ {xl , ac2 , at3 } , xX 

x3 = ac2 / . {0 -> 1. 1 -> -IJ] 
Xn[e3J := fx-uneEncoded = GSMDl££Eaicodeax*rak]ne [fz^me] 



> existing symbol Trajne" . 



Ouc[B3J~ {-X. 1, 



1, X, 1. 1, 1, 1, 1. -1, ■ 
, 1, -1., -1. -1, 1, 1, -1 
1, 1, -1. -1, -1, 1, -1, 



1. -1, -1, -1. -1, 1, 1. -1, -1, 1, _1, _1, _i, 1, 

1. 1. 1, -1. -1, 1, 1, -1. -1, 1, -X, -1, -1. 
-1. -1. -1- 1. -1. -1, 1, -1, -1. -1, 1, -1. 1. 1, 



xn[Bai -= MoaOutputrr-ame = HodulBtoz- [Z.] [ CDMAEncode [ f rameEncodod, QeldaBqlxat // taBtl 



In[B6] = 



SafTO [-HodOutpixtPxTEujiiaSiicodedOSHIjlka 
« MadOut.putFx-ajne.m; 
Lengrtli [HodOutsiutiFz-aine] 

73440 



a-, HodOutpufcFrai 



InC87J:= TedtSat^GL ■ 



TcUco£HodOut:prutFzrBme, {50, 7300}] j 

: new symbol name 'TestData" is 



InfBBj.-'- OlosCData// {Re [#] , I«n[#])t // Txanepoao // 

t-iatPlot [#, PlotiJodLned -> Tm&, X>lokRajca?« -> Ail, AapE>c:lLllat;a.o -> 1]&; 




In the PrimitivcCDMA 



we need to specify the sample. In the CDN4A synchi 



we discover the sample. 
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[B9J.-= CllB«ACo«r««By»cliraiiiBer-[jao«30ut:smt_, GoldSeq:_. Tlire^liol«5L_, Ov^r-Sanmliner 1 -- 
Module C(arl, x3, x3, x4. xSState, " 
acECouzit:, xGMaacCoarx', sbq, x7ljpdatia, sc8, acS, xlO, xlX, xl2, acl3>, 
atl = MoaOutpuc; 
x2 = aoXdSeg/. 0 -> -X; 

x3 = Talsle[ •^ModCi, 2] , {i, 0, tiMfftll [acZ] - 3.) ] ; 

x4 = roldJ:.d.st:[PX\iB, 0, X2] // Rest //I-#t//at3#6; 
xSState =xa]Ee[xl, ( X>en5rrli[x2] ■» 1} Overs Eunpllag] ; 
xCComxt: = 1; 

ee«i = r»rc9p[xl, Ov-ezrSaiiiDXlnv ( X.ea9t:h [ac3 ] * X> ] ; 
x7DEMla.ce := Module [{). 
x8 = V&xrlClDoL[x5SbBCe, OveirSa]qpl±sg] // Tz-azispoee; 
xS = MaE>[{Drop[«, -X] . x4, Drop[«, 1] . xa)&. x8]^ 
xlO = Mapti Al>B[Im[»[tlJ]]]. 

Aba[RB[#[[l]]J] , Ay>H[Rer#[:[2]]]] , AbB£Im[#[[2]]]] }&, x9]; 
xll =l£(Max[xXO] > Tlureshold, T»rrow[{xXO, xficotmt, Tjr«©>], 

{xecouat:, xlO, xSMoxCorx- = SUuc [x6HaxCorr, xlO] , False}]; 
xecount: = xfiCouat * X; 

xsscace = Jo±a[Dxoptx5State, OverSaiiipliiisr] , Ta)ce[Beq. OvearsaaipXliigj ] ; 
sect = Dxrop [seot, Over'aampllng] ; 
xlX ]; 

xX2 = cat:ch [Table [xVupdatis, {i, X, I.oi»grtb[se<2] / ovez-Samolln? 
xX3 = If [IiBBC [xXa] = = = True, 

Cl3KUVFlAeSyiicbrosLiaBr-[HodOut:put, xl2[[2]], xd, C>v«rEBiopllns7] , 
{"Failed to Ceaxse SyzicluroxilBe" , False}]] 

XnOOj := COHAX'ijaesyzicluraxii. eer [UoaouCputi_r 

TlEuc-esboldCorrelatLloCouxifc^ Coxx-eXBt:i.ngSeix_, OTr©rSomplliiS7_] : = 
Module[{xl< x2, x3, x4, x5, x6, eeQ, x7lipdiatie, xeFl»e. 

xBSecond, xSFlral:, x9Seoond< xlOFisrsc. xXOSecond, xlX. xX2, xX3, xi«}, 
xl = X£ [TlureBboldCo3rrelatJ.oCouzit > 3, 

Ttix-esliolclCoz7-elat:i.oCouxic -3, TbareataoldCoz-z-eXatrloCouxii^] ; 
x2 = CoxrelaeijagSeiz; 

x3 = mrop £HodOut:put., xX Oveirsamplliig] ; 
xA = C£(MAFOslCloi>Fliidez-[x3, x2, OverSainpllJiigr} ; 
xS = CDUAJE>OBlclozurlxiaejr[Dz7c>p [x3, Z.BXii?t:li [x2] OvexrSanipllzig^] , x2, Oire3rSaittpXi.ri9] j 

x6 = CDUAPosltioxaFlxide>:[m:op[x3, 2 X.ex>gi:h[x2] cy^raz-Saniipliniff] , x2 , OvexSaii«>Xia?1 i 
x7 =. »OB±tloaAvBr-aaor-[{x4 t ] . ^SC [1] ] r «K[ [1] ] }] J 
xBFl»t = Drop[x3, (x7[[l]] - X) Oversaiqpllzi^] // 

PartlcAon[», OveirSBiniplln?] & // Transpose// # [ [x7 [ [2] ] ] ] 
xBSeoood = Drep[x3, x7t[XJ] Ov^x-Sanpllnff] // 

Partltioxi[#, Ovez'SainplJ.ziff] & / / Vx-anapose // « [ £x7 [ [2] J ] ] sy 
scSFixat = xBFirst: // Fa3nilt:loci[«, X.ezig^tti [x2 ] ] & j 
xlOFi-x-at: = Hap[FUBctloa[x, x . x2] , xSFlrat] ; 

xSaecond = xSSecond // Pa3rtiit:loii[#, l.Bngtita,[x2] ] i^y 
xXOSeeond = Hap [Function Ix. x . x2] , xSSeeond] ; 
DoModuXatojrr{xlOFlret, xXOSecond}] ] 

lnl31]:= CI3IIAFosi.ti.oxaFlnder[HodOut:pub_, CozTz-aletlngSeO— # Owz-SBiiqpllnL£r_] : = 

Module [{xSSCaCe, x6Count, se^, x7Update, xS, x9, xlo, xXl, xl2, xX3}, 
xSState = TakelHodOutput:, ( Z.er>grfclL [Conelat Jji?Ee<z] ♦ 1) Ovez'Sanvllns] ; 
x6Cousit =1; 

eeq = Drop [ModOutput, Ovex-Samplln? ( Zjen9rb[Co3rz-elatxxiBSe(z] -i- 1} ] j 
x7crpdat:e := Module [{}, 
x8 = Part:d.t:lon.[x5Stat.e, Ovez-Suapllxx?] // TxaxkBTpoao; 

x9 = Blap[(IXi:op[8, -1] . Coxz-elBtln^Saiz. I3r-opt«, 1] . Coz-xelatlng-Seq} xB] ; 
xSCo'unt: = xGCount + 1; 

xSState = JodLn[Sirop[x55t:at:e, OvexSawpllzig^] . Talce [seq , OvexSamplijciLs] ] ; 
Beq = Dxc>p[Beq, Ov^axBampUncr] ; 
x9 ]; 

XlO = Table [x7T]pdatEi, {1, 1, 10}]; 

XlX = MaplnaBxed[Kax[{Abs [Re[#] ] , Abs [Xm [«]]}] xlO, {3}]; 
xl2 = Haplndexed [Apply [Plus, #]&, xXX, {2}]; 
xl3 = PoBltion[xl2, Maxtxxa]] ] 

We need to define a position averager. "Wc far npw just take the first element 

xnfSBJ-= Poal1:±(3iiAvsxager[PoBieloBX.lBt_] := Flirsc [Poslt:i.exiX.xat:] 
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0\xcll01l=: DEModulatOE- [{{-157. 056 + 0.651379 1. - 0 . 691379 - 157 . 056 X , 

157.056 - 0.691379 I. 0.691379 + 157.056 1, - 157 . 0S6 + 0 . 691379 I } . 
{-169.734 - 0.510132 1. 0 . 51 0132 - 169 . 73 4 1 . 169.734 + 0.510132 1, 
-0.510132 + 169. 73fl I, -169.734 - 0.510132 I}}] 

In{102]:= Toms = CimACCiaxdBSyxLclix-ozilaez' [taodOutpurPz-ana llTixao\1K, 250 4]&, 
Boldaealiat^ // Last , 100, 4] 

aut{102l= DeKo(auls.torH{l-02745 + 15G.6B2 I. -157.056 + 0.691379 1, - 0 . 691379 - 157 . 056 I , 
157.056-0.691379 1, 0.69137 9+157.056 1, -157,056+0.691379 1, 

-0.691379-157.056 1, 157.25-0.929033 1, 19-1272-6-93509 1, -7.37139-18 6305 1 
-18.8305 + 7.37139 1, 7.37139 + 18.8305 1. IB . 8305 - 7 . 37139 I , -7-57128-19.05181. 
-1.02745 - 15 6.6B2I, 157.25 - 0. 32 9033 I, 19 . 1272 - 6 . 93 509 I , -7 . 57 12 B - 19 . 05 IB I . 
-1.2651 - 156.875 1, - 6 . 935 09 - 19 - 1272 I . - IB . 83 05 + 7 . 3713 9 1 , 7 . 57128 + 19 . OSIB I. 
1. 02745 + 156.682 X, -157 . 25 + 0 . 929033 X, -19.1272 + 6.93509 1, 7.57128 + 19.0518 1 
1.2651 + 156. B75 I, 6.93509 + 19.1272 1, IB . 33 05 - 7 . 37 139 X . -7.57128 - 19.0510 1 
-1.02745 - 156.682 1, 157 . 25 - 0 . 92903 3 1 , 19.3484 -7.13498 1. 156.682 - 1.02745 1, 
0.691379+157.056 1, -157.25+0-929033 1. -19,1272+6.935091, 7 . 5712 8 + 19 . 051 8 X 
1.02745 + 156.682 1. -157.25 + 0.929033 1, -19,1272 + 6.93509 1. 7 . 57128 + 19 . 051B I. 
1.2651 + 156.875 1, 6.93509 + 19.1272 1. 18.6305 - 7.37139 1. -7.57128 - 19.0518 1 
-1.02745 - 156.682 1, 157.25 - 0.929033 1, 19.3484- 7.13496 1. 156.875 - 1.2651 1 
19.1272 - 6.93509 I. -7.37139 - 18.8305 1, -19.0518 + 7.57128 1, - 156 . B75 + 1 . 2 651 X , 
-19.1272 + 6.93509 1, 7 . 3713 9 + 18 . 8305 1 , 18.8305 - 7.37139 1, - 7 . 57128 - 19 . 05 IB I , 
-1.2651 - 156.875 X. - 6 . 93 509 - 19 . 1272 I , -19.0518 +7.57128 1, 
-156.675 + 1.2651 1, -19.1272 + 6.93509 1, 7 . 37139 + 18 . 83 05 I , 

19.0516 - 7.57128 1, 156 . B75 - 1 . 2651 1 . 19 . 3 4B4 - 7 . 13498 1 , 156 . 682 - 1 . 027 45 I , 
0.591379 + 157.056 1, -157.056 + 0-691379 1}. {-0.476141 + 169.7151, 
-169-734 - 0.510132 I, 0 . 510132 - 169 . 734 1 . 1 69 . 73 4 + 0 . 51 0132 I. 
-0.S10132 + 169-734 I, -169.734- 0.510132 1, 0.510132 - 169.734 1. 

169. B13 - 0.714545 I, 20.5326 - 6.8667 1, - 6 - 8 871 - 20 . 502 6 I . - 20 - 502 6 + 6 . 8 871 I , 
6 . 8871 + 20 .5026 I, 20 . 5026 - 6 . 8871 I . -7.42576- 21.3328 1, 0.479141-169.715 1. 
169. H13 - 0.714545 X, 20.5326 - 6.8667 1, -7.42576 - 21.3328 1, - 0 . 746535 - 169 . 794 I 
-6.6667 - 20.5326 I. -20.502 6 + 6.8871 I, 7,4257 6 + 21.3328 I. -0.47 8141 + 169.715 I, 
-163.613 + 0.714545 1. -20.5326 + 6.8667 1, 7 . 42576 + 21 . 3328 1 , 0.746535 + 169.7 94 1, 

6- 6667 + 20-5326 1, 20 . 5 026 - 6 . B871 I . -7 . 42576 - 2 1 . 3328 1 . 0.478141- 169.715 1, 
169.S13 - 0.714545 1. 21.3629 - 7.40536 1, 169 . 715 + 0 . 478141 X , 

-0.510132 + 169.734 X, -169 - B13 + 0 . 714545 I, -20 . 5326 + 6 . B667 I , 

7- 42576 + 21.3328 1, -0-476141 + 169.715 1, -169.613 + 0.714545 1, 

-20-5326 + 6. B6S7 I, 7.42576 + 21.3328 1, 0 . 746535 +. 169 . 794 I, 6.0667 + 20.5326 1, 
20.5026 - 6.8871 r, -7.42576 - 21.33 28 1, 0 . 47 8141 - 169 . 715 1 , 169-613 - 0.714545 1, 
21.3629 - 7.40536 1, 169.794 - 0.746535 1, 20 . 5326 - 6 . 8667 I, -6.8871 - 20.5026 1 
-21.3328 + 7.42576 I, - 169 . 794 + 0 . 746535 X , - 20 . 5326 + 6 . 8667 1 , 6-8871 +20.5026 1, 
20. 5026 - 6. 6871 X. -7.42576 - 21.3328 1, -0.746535 - 169.794 1, - 6 . 8667 - 2 0 . 532 6 X . 
-21.3325 + 7.42576 1. -169.794 + 0.746535 1. - 20 . 5326 + 6 . 8667 I , 

6.BB71 + 20-5026 I, 21 . 3 32B - 7 . 42576 1 , 169.794 - 0.746535 1. 2 1 . 3 629 - 7 . 4053 6 X , 
169-715 + 0.478141 1, - 0 - 510 132 + 169 . 734 I , -169.734- 0.510132 1}}] 

ltitX05i:= Seal = Table [I * (1) (-X)''(l + 1). {1, 1, I^oxxgtli [TomS C [1, l]]])] 

DucriD5;= {I, 1, -X, -1, I, 1, -I, -1, I, 1, -I, -1, I, 1, -I, -1, X, 1, -I, -1, I, 1, -I, 
-1. I, 1, -I, -1. I, 1, -I. -1, I, 1. -I, -1, I. 1, -I, -1, I, 1, -I, -1, I, 
1. -I. -1, I. 1, -I, -1, I, 1. -I, -1, X, 1, -I, -1, I, 1, -I, -1, I, 1, -I, -1, X. 1} 
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■-—>^106J:^ Toilv5[[l, 1] ] Seal 

bucfJOffJ=- {-156. 6B2 + 1.02745 I, -157 . 056 -i- 0 . 691379 I , - 157 . 05G + 0 . 6913 79 1 , 
-157.056 + 0.691379 1, - 157 . 056 ♦ 0 . 691379 I , -157.056 + 0.691379 1, 

-157.056 + 0.691379 1, -157.25 + 0.929033 1, 6.93509 + 19.1272 1, -7.37139 - 18.6305 1 
7.37139 + 18.8305 1, -7.^37139 - 18-8305 1, 7 . 37139 + IB . B3 05 I, -7 . 5712B - 19 . 0518 1 
-156.682+1.02745 1, -157.25+0.929033 1, 6.93509+19.12721. -7.57126-19.05181, 
-156. B75 + 1.2651 I, 6.93509 + 19.1272 1, -7.37139 - 16.6305 1, 7 . 5712 B + 19 . 051B I , 
156. 6B2 - 1.02705 I, 157.25 - 0.929033 1, -6.93509 - 19.1272 1. 7.57128 + 19 0516 1 
156.875 - 1-26511, -6.93509 - 19.1272 1, 7.37139 + 18.9305 1, -7.57128 - 19.0519 1, 
-156.682 + 1.02745 I, -157.25 + 0.929033 1. 7.13*196 + 19.3484 1. 156.682 - 1.02745 1, 
157.056 - 0.691379 1, 157.25 - 0.929033 1, -6.93509 - 19.1272 1. 7.57128 + 19.0518 1, 
156.682 - 1.02745 1. 157.25 - 0.929033 1, -6 . 93509 - 19 . 1272 I, 7.57128 + 19.0516 1 
156.675 - 1.2651 1. -6.93509 - 19.1272 1, 7.37139 + 18.8305 1, -7.57128 - 19.0518 1, 
-156. 6B2 + 1.02745 I. -157.25 + 0.929033 1, 7.13498 + 19.3464 1, 156 . 875 - 1 2651 1 
-6.93509- 19.1272 1, 7.37139 + 19.8305 1, -7.57128 - 19.0516 1. -156.875 + 1.2651 1 
6.93509+19.1272 1, -7.37139-18.83051, 7.37139+18.83051, -7.57126-19 05161 
-156.875 + 1.2651 X, 6.93509 + 19.1272 1, -7.5,7128 - 19.0516 1, - 156 . 875 + 1 . 265 1 I . 
6.93509 + 19.1272 1. -7 .37139 - 18 . B3D5 I, 7.57128 + 19.0518 1, 156.875 -1.2651 1, 
-7 .13498 - 19.34B4 I, -156.682 + 1.02745 1, -157.056 + 0.691379 1. 
-157.056 + 0.691379 1) 



InflJ4J.-= dsmodLfzTBin^ = Sign [Re [ O^oaiiS [ [1, IJ] SaqZ] ] /. {-1 -> X, 

Out[114J= {1, 1. 1, 1, 1, 1. 1, 1, 0. 1. 0, 1, 0. 1, 1, 1. 0, 1, 1. C 

0, 1, 0, 1, 1, 1. 0, 0, 0, 0, 1. 0, 0. 0. 1, 0. 0, 1, 0, 1 

1, 1, 0, 1, 0, 1. 1, 0, 1, 1, 0. 1, 0, 0. 1, 1. 1. 1} 



The first and last bits have been lost in the processing 



Out (1091= {1. 1. 


1. 1. 1, 1, 


1, 


1. 1. 


0, 


1, 


0, 


1, 


0, 


1. 


1, 


1, 


0, 


1. 


1. 0 


0, 0, 


1, 0, 1, 1, 


1, 


0. 0, 


0, 


0, 


1. 


0. 


0, 


0, 


1, 


0. 


0. 


1, 


0, 1 


0, 1, 


1, 0, 1, 0, 


1, 


1, 0, 


1, 


1, 


0. 


1, 


0, 


0, 


1, 


1. 


1, 


1, 


1} 



Out 1122}= 72 



inixX51:= trvmcatedf raioe = frame // Drop [», l]fc//Drop[», -1]6 

Cue I lis J = {1, 1, 1, 1, 1, 1. 1, 1, 0, 1, 0. 1. 0, 1, 1, 1, 0, 1, 1. 0, 1, 0; 0, 0, 1, 0, 

0, 1, 0, 1, 1, 1. 0, 0, 0, 0, 1, 0. 0, 0. 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 

1, 1, 0. 1, 0. 1, 1, 0, 1, 1, 0. 1. 0. 0, 1, 1, 1, 1} 



in[116J := demoilfxrens - truncRtedframe 

Out 1116}= {0, 0, 0. 0, 0, 0. 0, 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0, 0, 0, 0. 0, 0, 0. 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0. 0, 0, 0. 0, 0, 0. 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0, 0, 0, 0, 0, 0. 0} 



■ References 

. 1 'fiinary Se^OMOcmm wltb aold-I.i.Ice Conco2.».t.i.azx but Zisjrg'sxr Z.lnear' Spaja" "kjy 
Eez-daz- Bozt AS and P . Vl jay Kiimaxr XESE Tzrajia . on Ixi£oinna.tilaiitPlii90xy Vol 40 No 2 , 
Baarclil9B4 
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Rln^lvlBloat fusa^e = -RlngDivdalontHoduloRlsg] [ {r>*&c»ii±ii&cor,lfuxn«ra.t:or) ] = 

Cauot:lt:e&t:,B.dmazxi.d£iir,D«nomdJiAl:e>sr, Fla^9> x£ xuun ^ {X, 0 , 0, 0, 0 , 0 , 0. 0, 0 , 1, 

0, 0,0, 0, 0,1} , den = {1,0,0,1.0,1}, cben alzi0Dlvi.8lon[2} [ {d^,ziuin} ] = ({1, 
0. 0,1, 0, 1, 1, O, O. 0,1}, {1, 0,1, O, 0} , {1,0, a,l,0, l>,rriaio} . tAben tlx© 

algoxltilun encount^eiTB ^ci z«ro dovlsor tlie flag- 1* Co Palea. a.sr RlasrSi.-ir±elozi[ 

4] [{{2,1,1,1}, {1,0. 0,1}}] = {{}. {1. 0,0,1}, {2, 1,1.1}. Falae} " 

Dz-opI.eadi^grZex'os : : visage = 
"I>roi;>I.eadlngZeroB [DenorolnatogSeg] Oxoeib tlie leaOlns zeros In t:be se^axxca, 
o.ff r»ropI.»a«aixiHZBroii[{3,3,l}] = {3,3,1}. DropLoacHagZwx-oB [ {0, 0, 0} ] o {}- 

PoadiJsleDlvlsazrd: zuaag^^ = "»aaailileDivia<ir<a [Modultfillljag] [Ox-dex-] givBd the list o£ 
possible devdLeoirs. Kot:^ t:he MS pcnraz' o£ Cbe poly ia UI5. Pi:>asililQZ>i-vJ.socB [4] [ 
2] = {{1},{2},{3}. {l,0},{l.X},{l,2},(l,3},{2,a},{2,l},{2,2},{2,3},{3,0}, 
{3,1}. {3.2}, {3.3}, {1,0,0}, {1,0,1}. {1.0, 2}, {1.0. 3}, {1, 1, 0} . {1, 1, 1} . {1, 
1,2}, {1,1, 3}, {1,2,0}, {1,2,1}, {1,2,2}, {1,2,3}, {1,3,0}, {1,3, !},{!, 3, 2}, 
{1.3,3}, {2,0.0}, {2,0,1}, {2.0,2}, {2 , 0, 3} , {2, 1. 0} . {2 . 1. 1} , {2 . 1, 2} . {2,1, 
3}, {2,2.0}, {2,2,1} . {2,2,2}, {2.2,3}. {2,3, 0}, {2.3,1}. {2,3,2}, {2.3,3}, 
{3.0,0}, {3,0,1}, {3,0,2}, {3.0,3}, {3,1,0}, {3,1,1}, {3,1,2}, {3,1, 
3}. {3, 2.0}, (3. 2,1}, {3, 2. 2}. {3. 2. 3}, {3. 3.0}. {3, 3.1}. {3, 3, 2}, {3.3,3} }- 

PolynomlaUHult:l.plicatilcmz susafre = "PolynomlalHult:lpllea1i3.on[IBeauloIl±ag] [ 
polyl.poly2] p«r-£oxma polyziomial mulirlpllcadoxi aoduloRi.119 m.g. 
Volynoiii±alBiult±plicotion[4] [{2.1}, {2,1}] = {!}- 

ModuJ-oMultpllcatioxi: lUSag^Q = 

■■UoduloICulfcpll.catloxxfModuloai.zi9} [Px-lndLtivePolyiioinial] [polyl,poly2] o. 
3 Kodulo«ulfcpXicatloB[4] {{1,2,1,3}] [{3.3,0}, {3,3,0}} = {1,0}- 

Rj^gPoHcz-r ru^aga = "lU^SPcnirex' [RljagHoduloJ [PxrliB±ClvsKl«iDenc] {BleiiMaiC] ZMnerPotrar^ [ 

4] [{1*0,0,3,2.3}} [{1,0,0}] = {{1.0,0}, {1,0, 0,0,0}, {1,2, 1,0}, {2, 1,1, 
2,1}, {1,0,2,0,1}, (2,1,3,1,0}, {3,3,1, 0,1}, {1,3,2.1,3}. {2,2,0,3, 
3}, {1,1.2,2}, {1,2. 3, 2,1}, (3,3,1,1,2}, {1,0,3,1,3}, {3,2.1,1, 
0}, {1. 0. 0.3,2}. {1,0}, {1,0.0,0}, {1.2. X}, {1.2,1,0,0}, {1,1. D, 
1,2}, {2,1,3,1}, {1,3.3,0,2}, {3,1,3.3,3}, {3, 2, 2. 1. 1} , {2. 0, 1, 
3,2}, {1,1, 2, 2.0}. (2, 3, 3, 3.1}, {3.1.0,0,3}, {3.2,1,1}, {2,1, 0,2,3}, (1}} - 

CyolioMttltJLplati-wia(?raup: suaago =."CyellcMiiltiplatlv©Group[ModuloRlng-] [ 

Pz-lulkl-vePolyxiaiidLal} e.g. Cy«:la.cHiilCi.pl&t:i.vearoup[4] [ {1. 2, 1, 3} ] □ {{1}.{1, 
O}, {1,0,0}, {2, 3,1}, {3,3,2}, {1,3,3}, {1,2,1}} and Cycllc!iaal1:iplati-\reOx-oup[ 
4} [{1,0,0,3,2,3}] = {(!}, {1,0,0}, {1,0, 0,0,0}, (1,2, 1.0}, {2,1,1,2,1}, 
{1,0,2, 0,1}, {2,1,3,1,0}, {3.3,1, 0,1}, {1,3,2,1,3}, {2,2, 0,3,3}, {1,1,2, 
2}, {1,2,3,2,1} , (3,3. 1.1,2}, {1,0,3,1, 3}, {3,2,1, I. 0}, {1,0,0,3,2}, 
{1,0}, {1,0, 0,0}, {1,2,1}, {1,2,1,0,0}, {1,1,0,1,2}. (2, 1.3,1} - 
{1,3. 3, 0.2}, {3, 1.3, 3. 3}, {3,2.2.1.1}. {2,0,1.3.2}, {1.1.2,2, 
0}, {2, 3, 3. 3, 1}, {3, 1, 0, 0, 3}, {3,2,1.1}, (2, 1,0. 2.3}} - 

zeroSsqruonceat ;uBa?e - 
"Zex-oBeciuencea [ModuloKln?] [Order-] genez-ated tta zez-a Ideal zar-ose^enaeia [4} [ 
3]= { {0.0,0}, {O, 0,2}, {0,2,0}, {0,2, 2}, {2, 0,0}, {2. 0,2}, {2.2. 0}, {2,2,2}} - 

ZeroPad: tusag-e = -ZeroPad[Oi:dex-] [zlenttenti] ejf zexoPad[3] [{1}] = {0,0,1}- 
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, ) T-»oAa±cETO«llBioa[4] [(1,2,X,3}] = {{(0,0,O),{D,O,O},{}}. {{0,0,0), to. 0,1}, 

{2>>,{{0.0,0>. {0,1,0}. {2,0)}.{{0.0,0>,{l,0,0>,{2,0,0>},{{O,O,oi;{2.3, 
X>,{2,2>).f {0,0,0}. {3.3.2}. {2,2.0>>,{{0.0,0},{X.3,3).{2,2.2)}, {{0,0.05 
{1.2.1}, {2, 0.2}}, {{0,0,1}, {0.0.0}, {!}},{ (0.0.1}. { O . 0 , 1} , {3 } } , { { 0 , 0 , 1 > 

1},{3,3.2}, {2,2.1}}. {{0,0,1},{1,3,3},{2.2,3}},{{D,0,1}, (1,2,X},{2,0, 
3}},f {0.1,0}, {0,0,0}, {1,0}}. ({0.1,0}, {O, 0,1}, {1,2}}, {{0.1,0}, (0,1,0}. 
{3»0}},{{0,1,0}. {1,0,0}.{2.1.0}},{{0,1,0},{2.3,1},{3,2}},{{0,1,0}.{3.3, 
2}, {2, 3,0}}, ({0.1,0}, {1.3, 3}, {2, 3. 2}}, { { 0. 1, 0} , { 1. 2, 1} , (2, 1. 2} } . { {1, O , 
0}. {O.O.D}.{1, 0,0}}, {{1,0.0}, {O, 0.1}, {1.0, 2} },{{!, 0,0}, {0,1,0}. {1. 2' 
0}}. ({1.0,0}, (1,0,0}, {3,0,0}}. {{1,0,0}, {2,3,1}, {1, 2, 2} } . { (1, O, 0} , (3 , 3 , 
2}, {3,2.0}}. {{a*0,0},(l,3,3},{3,2,2}}. { {1,0.0}. {1.2,1}, (3, 0,2) }.{(2l 3' 
1},{0,0.0}.{2,3,1}}, {{2,3,1},{0.0,1},{2,3,3)},((2,3,1},{0,1,0},(2,T, 
1}}.{{2,3,T},(1,0,0},{3,1)}. {{2.3,1}, (2, 3, X}, {2 , 1. 3 } } . { (2 , 3 , 1} , (3 , 3 , 2 } , 
(1,1}}, {{2. 3,1), (1,3, 3}, {1.3}}, {(2. 3.1), {1.2,1}, {3. 3)}. {{3, 3, 25, {0.0, 
0}, {3,3,2}}. {{3.3,2}.{0,O,l},(3,3,0}}.{(3,3.2},{O.l,0},(3,l,i}}:{{3,; 
2}, {1.0,0}. {1,3, 2}}, {(3, 3, 2}, (2. 3.1}. {3, 1.0}), {(3, 3, 2). {3, 3. 2}; {l,i: 
2}}. {(3,3,2), {1,3, 3), {1,1.0}}, {{3, 3, 2}, (1,2,1), (1,3.0}}, ({1,3,3}, {0.0, 
a),{l,3.3)},{{l,3.3},{0,D,l},(l,3,l}),{{l,3.3}.{0,l,0),{l.l,3}},{(l 3! 
3). (1.0,0), (3, 3, 3}), {{1,3, 3}, (2, 3.1), (1,1,1)}, {(1,3. 3), (3, 3. 2}, (3,1; 
3}}.{{1,3,3},(1,3,3},{3.1.1}}. ({1,3,3}.{1.2,1),(3,3.1}},{{1.2,1).(0^0. 

0) . {1.2.1}}. {{1,2,1}, (0,0,1), (1,2. 3} },{{1. 2, 1},{0,1,0},{1.0.1)}: {{1^ 
2, 1}, (1.0,0), {3, 2,1)}, ({1,2,1}, (2. 3,1). {1,0, 3)}, {{1,2, 1), (3, 

3,2}, {3. 0.1}), ({1.2.1). (1.3. 3), {3,0. 3}}, {(1,2.1), (1, 2, 1} , (3 , 2. 3} } } - 

ozsuBsg-a = "o [MoauloBJLnff] [Prlinltlv-Molynonilal] [TwoadLLcEXeonoiit] •-g - 

Atttomor phH rnnn lgmar rtiaaga 'Alit:aanoxp]alBmSlgiii& [HoduloRijag^} [Pjrlmlt:lv«Valyiuimxal] • a 
AiiteBiiorplilsmHlgiiia[4] [(1.2,1,3)] = { ( {0. 0, 0} , (D, 0, 0) } , ( (0, 0,2), (0, 0. 2) ) , 
({0,2,0}. (2,0,0)), {(2,0.0). {2.2,0}}, { { □. 2, 2) . (2, 0, 2} } . { (2. 2. 0} . (0 , 2, 
0}}, {{2, 2. 2}. {0,2, 2}), {(2,0, 2>, {2. 2 . 2} } . { (0. 0. 1) . { 0 , 0, 1} } , { (o, 0, 3) . (0 
0,3}}. {{0.2,1}. {2,0, 1}},{ (2, 0,1). (2, 2,1}}. ((0.2, 3). {2,0, 3)), ({2,2,1}. 
{0,2,1}), ((2, 2, 3), (0,2, 3}}, {{2, 0,3), (2, 2, 3}}, (CO, 1,0). {1,0,0}}. ({0,1, 
2},{1,0,2)}, {{0,3,0}, (3. 0,0}}, ((2, 1,0}, (3, 2,0}), ({0.3,2}. (3, 0.2}}, ((2, 
3.0}, {1,2,0}). {(2. 3, 2), {1,2, 2)), ({2, 1,2), (3,2,2)), {{1, 0, C) , (3, 3. 2) } , 
{(1.0,2), {3,3,0}}, ( (1,2,0), (1,3,2)}, {{3,0,0}, (1,1,2}}, {{1,2,2}, (1,3. 
0}). {{3,2,0). {3,1,2}), ({3,2.2}. {3.1,0}}. {{3,0,2). (1.1.0)}, ({2,3,1}. {1, 
2,1}), ((2, 3, 3), (1,2, 3}}, {{2.1.1}. {3, 2,1}}, ({0,3,1}, (3, 0,1}}. {(2,1,3), 
{3,2,3}). ({0.1,1}, {1,0,1}}, { (0,1,3), (1,0,3)}. ((0.3.3). {3,0,3}}, {{3,3, 
2}, (0,1,0)}. {{3, 3,0}, (0,1, 2}}. {{3,1,2}. {2.1.0}}, ((1,3, 2), (2, 
3r0)},( (3, 1,0}. (2, 1,2)}, ({1.1. 2), {0,3,0} ),{(!, 1,0), (0,3, 2}}. 
{{1,3,0), {2, 3, 2}}, {(1,3, 3}, (2, 3,1)}, ((1,3,1), {2, 3. 3}}, {{1,1, 
3). {0,3,1}}, {(3, 3, 3), (0,1.1)). ({1.1.1}, (0.3, 3)}, {(3,1,3), (2, 
1,1)}, {{3, 1.1), (2. 1,3}}. {(3,3,1), {0.1,3}}.{{1.2.X}. (1,3,3)). 
{{1,2,3}.{1,3.1}},{{1,0,1},{3,3,3}},{(3.2,1}.(3.1.3}},{{1.0, 
3}, (3,3,1)) . {(3,0,1), (1,1,3)}, {(3,0,3), {1,1, 1)}. {(3,2,3), (3,1,X})) - 

UnlksKlng; : usage = "iniltisRlxig [HoaxLloRi-ns] [E>riiid.^ivaFolyziaiBd.sl] e.g. njoitsKlxigr^l [ 
(1.2.1,3)] = {(0,0,1), {0,0,3). {0,1,0), (0,1.1), {0,1,2}, (0,1,3), (0,2, 

1) , {0.2.3}. {0,3. 0), (0,3.1}, {0,3,2}, (0,3,3), (1,0, 0), {1,0.1}, {1,0,2}, {1, 
0.3). {1.1,0} ; (1,1,1). (1,1,2). (1,1,3), (1,2,0), (1,2,1), (1,2,2), {1,2,3}, 
{1,3.0), (1,3,1), (1,3,2), (1,3.3), (2.0.1). {2,0.3), (2. 1,0), (2, 1,1), (2,1, 

2) , {2,1.3}. {2,2.1}, {2,2,3}, (2,3,0), (2, 3 , 1) , (2, 3, 2} . (2. 3 , 3 ) , ( 3, 
0,0), (3,0,1), {3.0.2}, (3,0,3), (3,1,0), (3.1,1). (3,1,2). {3.1, 

3) . (3. 2.0}, (3,2,1), (3.2,2). (3.2,3). {3.3.0}. (3.3.1). (3,3,2), (3,3,3} }- 

r-zyx0ago -7-[Oxd€>r_] [er_] [»:_] e -gT-E*] [x3] [{3, 0, 3)] = - 
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E-srlniitivoPolynomialJ 0.3. !£r-ac«Repx«««ncatlon[4] f {1, 2, 1, 3 } ] = f f {0, 0, 0) , {0, 
0.0)}, {{0,0,2). {0,0,0} >, {{0,2. 0), {2. 2,0)}, {{2,0. 0>, {0,2,0)), {{0,2, 2}, 
{2,2,0}), {{2,2,0}, {2,0,0}}, {{2,2, 2>, {2,0,0}), {{2,0,2}, {0,3,0}), {{0,0, 
T}, {0,0,2}}, {{0,0,3}, {0,0,2}), {{0,2,1), {2,2,2}}, {{2,0,1}, {0,2,2}}, {{O, 
2,3), {2,2.2}), {{2. 2. X), {2. 0.2}). {(2, 2, 3), {2, 0,2)), {{2.0,3}. {0,2,2}}. 
{{0,1. 0). {1.1.0} }. {{0.1.2}, {1, 1,0) ), {{0,3.0), {3. 3,0}}, {{2.1,0). {1,3, 

0) }. {{0,3. 2), {3,3.0} )V{ {2.3.0 ), {3, 1,0)), {{2,3,2}, {3, 1, 0} ) , {{2, 1, 2} , {1. 
3.0)}, {{1.0.0}. {0.3.2}), {{1,0,2), {0.3,2}}. {{1,2,0}, {2,1.2}}, {(3,0,0). 
{0.1,2}}. {{1.2.2}. {2. 1,2}}, {{3, 2.0). {2,3.2}}. {{3, 2, 2} -{2. 3, 2}). {{3,0. 

2) , {0,1,2}}. {{2,3.1}, {3, 1,2}}, {{2, 3. 3). {3.1.2}), {{2,1,1), {1,3,2}}, {{O, 
3,1}. (3, 3, 2}). {{2, 1,3), {1.3, 2}}, {{0,1,1), {1.1,2}}, {{0,1,3), {1,1,2}}, 
{{0.3,3}. {3,3,2}). ({3,3,2}. {3 , 0, 2) } , { { 3 , 3 , 0} , {3 , 0 . 2} } , ( {3 . 1, 2} , {1, 2 . 
2}}. {{1.3. 2). {3, 2. 2}}, {{3. 1,0). {1.2. 2)), {{1,1, 2), {1.0, 2}). {{1.1,0}, {1, 
0,2} },{{!, 3,0), {3, 2, 2)), {{1,3, 3), {3, 2,0)}, {{1,3,1}, {3, 2,0)), {{1.1, 

3) , {1,0,0}}, {{3,3.3), (3,0,0)}, ({1,1,1}, (1,0,0)}, {(3,1,3), {1, 
2.0}}. {{3. 1.1). {1.2,0}}. {(3.3.1), {3, 0,0)). ({1,2.1), {2, 1,0}), 
{(1,2,3), {2,1,0}), {{1,0,1}, {0,3,0}}, ({3, 2,1), {2, 3,0)), {{1,0. 

3}, {0,3,0}), {{3, 0,1), {0,1,0} }, {{3,0, 3), {0,1,0}), {{3, 2,3), {2,3.0} )}- 

!PupleReip3:asezit:at:iozi: :uaa£re = "rrupleRepx-asezicat:lozi[HodiiloRln0] [Pr-li&lelvaPolyxuMni.al1 
deaitr&ed fox fielOs. e .g-.TupleRopxes«nt«t:ioii[2] [ {1. 0. O. 1. 0, 1) ] = {{0.0.0. 
0.1),{0,0,0.1, 0), {0,0,1, 0.0). (0,1, 0,0,0}, {1, D,O,0,O), {O, 0,1, 0,1), 
(0,1,0,1, 0), (1.0,1, 0,0), (0,1,1.0,1), (1,1.0,1.0). (1.0,0,0,1), (0,0, 
1,1,1), {0,1,1,1,0), (1,1,1,0,0), {1,1,1,0.1). {1,1.1.1,1}, {1, 
1,0,1, 1), (1,0,0.1.1). (0.0,0,1,1), (0,0,1,1.0}, {0.1,1,0,0}, 
{1,1,0,0.0}, (1,0. 1,0,1), (0.1, 1.1,1), {1,1,1,1,0}. (1.1.0,0, 

1) , {1.0,1,1,1}, {0,1, 0,1,1), {1,0, 1,1,0}, (0,1, 0,0.1), (1,0, 0,1,0) }- 

Klzxlmunfoly! tuaa^a = -MiAlmumPely {HoduleKla?] (PxlnuLtlvePolyzioiidLBLl] [z-aotPowex]'[ 
-vctx] designBcl £ox fields. e. 9. UlxdLaiuiiiPaly [2 ] [{1, O , 0 , 1, 0 , 1} ] [3] [x] " 

GXB.e££eiietil3od.: iUBacra 3 -ara«££«Hat:2]iod[OF2BidLspolv] senexabea cl prdLmltlvQ 

polyaamlal in xun? z4 m.g Gxa«££eltotl3C>d[ (1, O , O, 1, o , 1} ] = (1,0,0,3,2,3) 



■ 5eiauexLcea8XLez-&t:ox [ModuloRlz>9] [Px^lsiJ.bl'wFoly] [IzxlbialCozid] {Se<lZ.Bxi^1:lx] 
o , g . Seilua»caGeiiQz:»t:axT[A] [{1, O, 0, 3, 2, 3)] [{0, 0, 0. 0, 2)] [31] = 
({2, 0, 0. D, 0. 2. 0. 0. 2. 0, 2. 2, 0, 

0, 2, 2, 2, 2, 2, 0, 0, 0, 2, 2, 0. 2, 2. 2, 0. 2, 0}, (0, 0, 0, 0, 2)} 

laltialconditioiisj juiiago = "laltialConcU-tioriB [Pirimlti'vePoly] e.^ Zxii.tilalCoxidi.t:loxis [ 
(1,0,0.3,2,3}) = {{0,0,0.0.1). (0,0, 0.0. 3), (0,0, 0,3,1}, {0,O, 0,2, 
3), (0,0, 2, 0,1), {0,0, 2, 0,3), (0,0,2,2.1), (0,0. 2. 2. 3}. (0,2.0. 
0,1), {0,2, 0,0, 3), (0,2, 0,2,1), (0,2, 0,2, 3), {0,2, 2, 0,1), (0,2, 
2,0,3) , {0,2,2,2, 1). {0,2,2,2,3}, (2, 0, 0, 0, 1} , (2 , 0 , O, O, 3 ) , (2 , 
0,0,2,l},{2,0,0,2,3},{2,q,2,0,l},{2,0,2,0,3}.{2.0,2.2.1). 
(2. 0.2. 2. 3), {2, 2, 0,0,1), {2, 2, 0,0, 3}, {2, 2, 0,2,1). (2. 2, 0.2. 
3), {2, 2, 2, 0.1), (2, 2. 2. 0,3). (2. 2. 2. 2.1). (2. 2. 2. 2, 3), (0,0, 0,0, 2)) - 

Oold£eaue&cet ;ui99g9 = "C3oldSBixuaiice [PxlxBlCi-v^ePoly] [ZzalbCondSunJaBX] e.g. 
Boiase^axxc»[{l, 0,0.3,2.3}] [33] = (1,0,0,0,0,1.0,0,1,0,1, 

1. 0. 0.1.1. 1.1.1, 0,0, 0,1,1, O. 1. 1, 1, 0, 1, 0) " 

Aut:acorx-ela£.lc3xiSegii«xica: tuaa^a = *'Aiit:ocorxBlati.oiiSeQcudxice[aoiasaauanca] a>g. 
Aut:c>ce>irx-ela.t.l.oiiSsizu8zic«[{l, 0, 0 , 0, 0, 1, 0, 0, 1, 0 , 1, 1, 0. 0 . 1, I. 1. 1. 1. 0. 
0.0,1,1,0,1,1,1,0,1,0)] = {31,-1,-1,-1.-1.-1,-1,-1,-1,-1, 
-1,-1,-1,-1,-1,-1,-1,-1,-1.-1,-1.-1,-1,-1,-1,-1,-1,-1,-1,-1,-1} - 

speclfledaolOsaq^axicesi tuaage = 'SpeclfiBcLGoldseguencas [Pxlmltiv«Poly_] [ 
Iluiid3exXjl.Bt:_] a.g. SpDci.£lBdGoldSeiiuezicas [ (1, 0. 0, 3 . 2. 3 } ] [ { 1, 2 , 3 , Irftst - 

1.1. aBt)] = {{0,0.0,0,0,0,0,0,0,1,0.0.0,1,0,1,1,0,1,1,1, 

0,0. 0,1,1, 1,0,0,1.1}, (1,0,0,0.0,1,0,0,1, 1.1,1,0,1,1,0, 

0,1, 0,1,1, 0,1, 1,1, 0,0. 1,0, 0,1), (0,1, 0,0,0, 0,1, 0,0,0, 0. 
1,1.1,0,0,0,1,0,0,1,0,0,1,0,1,0,1, 1,1,0), (1,1,1,1,1, 0, 
0,0,1,0,0,1.1, 0.0,0. 1.1. 1.1, 1. 0,1. 0,1, 0,1, 1.1, 1,1), {1, 
0,0,0,0,1, 0,0,1, 0,1.1,0,0,1. 1.1.1,1.0.0.0.1.1.0,1, 1,1, 0,1,0} }- 
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(1,0,0,3,2,3}] = {(0,0,0,0,0,0,0.0,0,1.0,0,0,1,0,1,1,0,1,1,1,0, 
0,0,1,1,1, 0,0,1,1), {1,0,0,0,0,1, Q, 0,1, 1,1, 1,0, 1,1, 0,0,1,0, 1,1, 
0,1,1,1,0,0,1, 0,D,1}, {0,1, 0,0, 0,0, 1,0, 0,0, 0,1, 1,1, 0,0,0,1,0,0, 
1,0,0,1,0,1,0,1,1,1,0}, {1,1,0,0,0,1,1,0,1,0,1.0,1,1,1,1,1,0,1, 
0, 1, 0, 1. 0, 0. 0, 1, 0. 1. 0^ 0} .{0.0,1, 0.0,1,0.1,1,1,0,1,1,0,1,0,1, O, 

l,o,o,o,i,i,o,i,o,o,i,i,i}, {1, o, 1, o, o, o, o.i, 0.1,1, 0,1, o,o,i,o, 

1,0,0,0,0,0,0,0,0,1,1,1,0,1}, {0,1,1,0,0,1,1,1.1,0,0.0,0,0,1,1, 
0,1, 0,1,0,0.1,0,1, 1.1,1,0.1,0), {1,1,1,0,0,0,1,1, 0,0,1,1,0,0,0, 

0. 1.0.1.1.0.0.0.1. 1.0. 0.0. 0.0.0}, {0,0, 0,1,0, 0,1, 0,1, 0,0, 0,1,0, 

1, a,o,a,i,i,o,i,o.l,o.a,i,i,o.D,i).{i,o,o,i, 0,1,1,0,0,0,1,1,1, 
0,0,1,1,1,0,1, 0,1,1, 0,0,1, 0,0, 0,1,1), {0,1, 0,1, 0,0,0,0,1,1,0,1, 
0,0, 1,1,1,1, 0, 0,0,1,0,0,1, 0, 0, 0. 1,0,0), (1,1, 0,1,0,1,0, 0,0,1,1, 
0,0. 0,0,0,0,0.1,0, 0,1,1,1, 1,1,1,1,1,1,0), (0,0.1,1,0,1,1,1,0, 0, 
0,1, 0,1, 0,1,0,0,1, 0,1,1,1, 0,1,0,0,1, 1, O.I), (1, 0,1, 1,0, 0,1. 1,1, 
0,1, 0,0,1,1,0,1,1, 0,0,1,1,0,1,1,1, 1,0,1,1,1), (0,1,1,1, 0,1,0,1, 
0.1, 0,0, 1.1, 0,0, 1,1, 0,1, 1.1.1,1, 0,0,1, 0,0,0,0), (1,1,1,1, 0,O, 0, 
1,1.1,1,1,1,1.1,1,0.0,1,1,1.1,0,0,0,1, 0,1,0,1,0), (0.0,0,0,1,0, 
0,1, 0, 0,1,0,0,0,1, 0,0,1,1,1,1,1,1, 0,0,0,0,0,1,1.0), (1. 0.0. 0.1, 
1,0.1,1,0,0,1,0,0,0.1,1,0,0.1,1,1,0,1,0,1,1,1,1,0,0), (0,1, 0,0, 
1,0,1,1,0,1,1,1,1,0,1,1,1,0,0,0,1,1,1,1,1,0,1.1.0,1,1). (1.1.0, 
0,1. 1. 1,1. 1,1, 0,0,1. 0.0,0, 0,1,1, 0,1, 1, 0,0,1, l,0,O,O, 0,1), (0,0, 
1,0,1, 1,0, 0.1, 0,1, 1,1, 1,0,1; 0,1. 1,0, 0,1, 0,1, 1,0, 1,0, 0,1,0), (1, 
0,1,0,1,0,0,0.0,0,0,0.1,1,1,0,1,0,0,0,0,1,1,0,1,1,0,1,0,0,0), 

(0,1,1,0,1. 1.1,0,1. 1,1,0,0.1, 0,0,1,0,0,1,0,1,0.0,0,0,0,1,1, 1, 

1), (1.1,1,0,1,0,1,0,0,1,0,1,0,1,1,1,0.1,1,1.0.1,1,1,0,1,1,0,1, 
0. 1), (0, 0,0,1, 1, 0,1, 1,1,1,1, 0,1,1,0,1,1,1,1,1.0,0. 1.1,1,1,0. 1, 
1,0.0), (1,0. 0,1,1,1,1,1,0,1,0,1,1,1,1,0,0,0,0,1,0,0,0,0,1,0,1, 
0,1,1,0), {0.1,0,1,1.0.0,1,1,0,1,1.0,1,0,0,0,0,0,0,0,0,1,0,0,1, 
1,0,0,0,1),(1,1, 0,1,1,1, 0,1, 0,0,0,0.0,1,1,1,1,1, 1,0,0,0,0,1,0. 
0,0,1.0,1, 1), (0,0, 1,1,1, 1,1,0,0,1,1,1,0.0.1.0,1,1,1,0.1. 
0,0,0,0,1,1,1.0,0,0), (1.0,1,1,1,0,1, 0,1, 1,0,0,0,0,0,1, 

0,0,0,0,1,0,1,1.0,0,0,0,0,1,0), (0,1,1,1,1,1,0,0,0,0,1, 

0,1,0,1,1,0,0,0,1,1,0,0,1,1,1,0,0,1.0,1), {1,1,1,1,1,0, 
0,0.1,0,0,1,1,0,0,0,1,1.1,1,1,0,1, 0,1,0,1,1,1.1. 1).(1. 

0,0,0,0,1, a,o, 1,0,1, 1,0, 0,1. 1,1. 1.1, 0,0, 0,1, i, o, i,'i, i,o,i,o}) - 

croseeerr»lacd.oziS«Queiic&; tuea^e = "Crosscorz-eladaziEagu«iie«[{ColdSeal,GolclSeq2) ] 
e. g^. aaqllBfc =Speci£ledGoldScH3UBiieea [{1.3,2.0,1, 1,1,1, 1, 1,1)]I{1, 2, 3, 
4,5,6,7, e,Z,a0t -l,X.a3t:) ] ;C>ut:er-[ (crosseoxx-elaclonSewe&ce [ (#1, #2) ] //Uxilaii ) u, 
Bea:li.Bt,fie(illet:,l] = ({(-33.-1.31,1023), (-65, -33 , -1, 31, S3) , (-65,-33,-1, 
31, 63), (-65, -33,-1,31,63), (-65,-33,-1,31,63), {-65, -33, -1, 31, 63 ), (-65, -33, 
-1,31,63), (-65,-33, -1,31,63), (-65, -33,-1,31,63) . (-33.31)). ((-65, -33,-1, 
31.63), (-33,-1,31,1023), (-65,-33, -1, 31, 63) , (-65, -33, -1,31, 63) . (-65,-33, 
-1,31.63), (-65,-33,-1,31,63). (-65,-33,-1.31,63), (-65,-33,-1,31,63), (-65, 
-33.-1, 31,63), (-33,31)), ((-65,-33,-1,31. 63) , { -65. -33, -1, 31, 63) , (-65,-33, 
-1,31,63,1023), (-33,-1,31), ( -65 , -33 , -1, 31, 63 ) , (- 65 . -33 . -1. 3 1, 63 ) , (-65, 
-33,-1,31,63), (-65,-33,-1,31, 63), (-65, -33,-1,31, 63) , (-33. 31) ) . {(-65, -33. 
-1.31. 63). (-65,-33,-1.31, 63). (-33,-1.31). (-65,-33,-1,31, 63,1023), (-65, 
-33,-1,31,63), (-63, -33,-1,31, 63), (-65,-33.-1.31, 63). (-65,-33,-1.31, 63), 
{-65,-33,-1,31,63), (-33,31)), {(-65,-33.-1,31,63), (-65.-33.-1.31.63), (-65. 
-33,-1,31,63), (-65,-33,-1,31,63), (-33 . -1. 31, 1023) , (-65 , -33 , -1, 31, 63) , 
{-65,-33,-1,31.63). (-65.-33,-1,31,63), (-65.-33.-1,31,63), (-33,31)), {(-65, 
-33,-1,31,63), (-65,-33,-1,31, 63), (-65. -33 , -1, 31, 63) , (-65. -33 , -1, 31, 63 } , 
(-65,-33,-1.31,63), (-33 , -1, 31, 1023) , (-65 , -33 , -1, 31, 63 ) , (-65, -33 , -1, 31, 
63). (-65,-33,-1,31, 63), (-33.31)), {(-65,-33,-1,31,63), (-65,-33,-1,31, 63), 
(-65,-33,-1.31, 63), (-65,-33,-1,31,63). (-65, -33, -1, 31, 63 ) , (-65, -33 . -1, 31, 
63), {-65,-33,-1,31,63,1023). (-33,-1,31), (-65,-33,-1,31. 63), {-33.31}), 
((-65.-33,-1,31.63). (-65,-33,-1,31,63), (-65, -33, -1,31, 63), (-65. -33,-1, 
31,63), (-65,-33,-1.31,63), (-65. -33,-1,31,63). (-33 , -1, 31) , (-65, -33, -1, 31, 
63,1023), (-65,-33,-1,31.63). {-33,31}), { (-65 , -33. -1, 31, 63) , (-65, -33,-1,31, 
63), {-65, -33, -1,31, 63), (-65,-33,-1,31.63), (-65,-33,-1,31,63), 
(-65,-33,-1,31.63), (-65,-33,-1,31, 63), {-65. -33, -1. 31, 63) , (-33, 
-1,31, 1023), (-33, 31)), {{-33, 31). (-33, 31), (-33, 31), (-33, 31), (-33, 
31), (-33. 31), (-33, 31), (-33, 31), (-33, 31), (-1,1023))) - 

Bagin. t " ' Px-lvnt • * " ] 
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ZlijagrZ>l.vlsi.oxi.[MoauloRlAg'_] [ {I3enami.iiatox'_, HuanBr'Bt:or_} ] 
Module [{xl, x2. x3, x4, x5, x6. Ox, Rx} , 
Cl<.ar[xl] I 

3c3.[i^. J_] := Mod [id, ModuloRing] ; 

x2 = a:ablo[atl[D»iioiiiiiiatox-[ [1] ) , J], {J, 0, ModluloRiiiff • 
x3 = CoCTpXemontL Table [i, {1, 0, WoauloRing - 1) ] , x2]j 



Module t{'uxlr ux2, ux3, uxA}, 

If [M«mb«r-Q[x3, Mx[[l]]], Tlixow[{Qx, Mx. Denewnlnacer-, ralao) ] ] 7 
uxl = TaJce [Rx, Ziengrt^h [Penomlnatoir] ] ; 

ux2 = Position [x2, 11X1 [[I]]] //Flatteja// (« [ [X] ] -1) ft; 
ux3 = noil[\Lxl — 11x2 PanomiTifttor, HodiiJ.ORi.ilB^] ; 

Qx = Jo±ii[Qx, {11X2}]; 
ixxA = Dx-op [NX, Zien£rt:li[nerkomlxiat:or] ] ; 
NX = aoijD. [Rest [ux3 ] , ux4 ] ] ; 
C«Cc1i[Tb1>1«[' No^lUpdate, {1, 1, I.eng«i [Numerator] - Length [Denomlimtor] + l}]j 
{Qx« Kx, IXTLOTn ■Lna t or< True}] 

] 

DropIjeadizLTZeiros [Danranl nn toi-Beq_] : = 
Module [{acl J, 
3H[{}] := {}/ 

»i[DS_] /r ual [1] ] ===0 j= xl[ Rest [DS] ] J 
xl [DS_] := I3S/ xl[ Deaomlnatoz'Seix] 
] 

PoBBll>leDdL-viBorB[Hod\iloSijag_] [02:d©3c_] : = 
Module [{xl, x2, x3, x4}. - 
xl = Table [ {x2 [±] , 0, ModuloHlng -1}, <1, 0, Ordew}]; 
x3 a Taljle [xa [i] , {i, 0, Order}]; 
x4 = APPly[Tal9le, Seizuence [ Jolii[{ac3} , atl] ]] // Fla.tt©xi[«, Order] & // Rest; 
Map[Drox>Z<eB.dlxi0Zeroe, »a] ] 

PolyiiaiiiialMultiplieatlon[ModuioRiag_] [polyl_, i>oly2_] : = 
Modulo [{3tl, »c2, 3€3, xA, atS, atS, x7, xB, x9}, 
xl = X.eia.9t2i[polyl] ■*■ Xian^b [poly2 ] - X; 
x2 = Tal>Xe[x3 "i, {1, O, xl -1>]; 

x4 = (Ta3co[3t2, I^ngtli [polyl] ] Reirerse) . polyl; 
xB = (Take[x2, I.exL9tlx[poly^2] ] //Reverse) . poly2j 

x7 = Table [Coefficiexit[x6, x3 " (xl - 1 -1)] // Kod[#, ModuloRijagJ ft, {±, 0, xl - ; 
xB = Joi.n[x7, {Mod[x6/. x3 -> 0. HoduloRln?] ) ] // urepLeadlngZezYia] 

We define polynomial multplication, modulo a primiiive polynomial 

Hoduloiailtplicatlon[HoduloRliig_] [PrimitiirePolyiiomial^] [polyl_, poly2_] : = 
Modulo t{xl, x2, x3, x4}, 

xl = PolynoiiilalKultipll.catlozi[ModuloRlsis] [x>o3.yl, poly2] ; 
x3 = Ri.nsrl3i.vl.8l,OBlHoduloRl2ag] [ {Prlmltl-vePolyiiomi.al, xl}] ; 
>^ [ [2] ] // mropXieadl&TZeros] 

RxxigFower [RlAgHodulo_] [Prlinit ^.v-eEleiuent^J [Element^] ; = 
Module [{xl, Bt, Noelupdate}, 
xl = i.eii£rth.[Prl2iiltiveElement] -1; 
Bt = {!},- 
Noelx^date : = 

St = Modnl oiynil tplloati.on[R±ngModulo] [Frlmltl-ireEleiBexit] [ at. Element] ; 
Table [Noelupdate, {±, 1, a**'-!}]] 
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CyclicMiilt:ii>a.ativ^«nip[ModiiJ.oRiiisf_J [PariiiLLti-vaPolyiiciniial_J s = 
Modu3.e[{x3., x2. x3, xA. xS. x6, x7, xB, x9 } , 



xl = (X«ng^la[PrlAltl-vaPoXync»nlal] ■•■ 1) / 2 / / Floox-; 
x2 = POHslblel)ivlBorfl[HodvaoRiiiB] [jci] ; 
x3 = taap[RiagrDlvlBleii[HeduleR±xig'] [ Px'lJiiJ.tlvePelyiioBilaJ.) ] 4t2]; 
= Select [ac3, »[r<]] =,= ,.*nio t] j 
atS = Select [x4, #[£3] ] == { 1> / / •TareuoepoBe / / aoia [# [ [ 1] ] , « [ [3 ] ] ] 

x6 = B«Xoct[x4, #[[2]] == {ModuioRiiasr - 1} &] // Tranapoae // Joia[» [ [1] ] , #[[3]]]6i; 
x7 = jroln[<{l}< {HoduloRlzka -!}}« x5< x6] //Union; 

xS = Map[ModuloMllltplleaclon[MoauloR±iia] [PsrimltlvePolysieamial] t». x?] // trjxion.; 

x9 = RijiePcFWBi:tMoAJ.oRixiB] CP>rJjna.tivePolync;iii«i»Xl [x8t[2] 3 ] // RofcateR^ ] 



ZQr'eBaq[ueaeaa[Meaulellljair_] [OsTdear..] : = 
UoduleCCxl, xa, x3, x4}, 
xl = Table ( {3c2 [1] , 0, 1>, <i, 1, Ordex- } ] ; 
x3 = TableE KoduloRizigr/ 2 x2 [1] , 1, Ox-der-}]; 

x4 = Ax>ply[*Bl>l»» Se<iuenca[aolaH5t3), xl] ] ] // Flatc»a[», Ordsr-llc ] 

ZoroPaa[Oiraer_] [Element^] = = 
Kodule[{xl, x2, x3}, 
xl = Table[0, {i., 1, Order - XiengchlElement:] } ] ; 
x2 = Joia[xl, El«B»nt] ] 

TMoAdJ.cEbKpBiii9d.oxifKi3duloRj.zig_] [Pz'JjnJ.td.-v-ePolyxiaiiiJ.al ] : = 

lSOdula[{xl, x2. x3, x4, x5, x6, x7, xB. x9}, 
xl = IjMi^h[Prliiilt:ivePolyxiosiilal] - li 

x2 = CV«=lJ.cMultJ.plBtJ.-TOG3rovLp[ModuLloRJjaB] [P2rij1u.tJ.-1rBP0ly11a11u.al] ; 

x3 = saap[ZeroPad[xl] . x2] ; 

x4 = Join [ {Table [0, {1, 1, xl}]}, x3]; 

x5 = rimctlonlx, ltoatx[ [1] ] ♦ 2 x[ [2] ] , ModviloRliiff] ] ; 
x6 = Table [{x« [[!]], x4[[j]], x5[{x4[[i]], x4[tD]])] // Dropl.s»adlnffZeroe} , 
<J., 1, I.anB*Hi:at4] >, ti . 1, I.enertbtx4] ) ] // 

FlactaxiCff, a]t] 

c [18od\ilollixi?_] [PrJjailt±-woPolyimmJ.al_] [T«roadi.CEleiBAnt:_] : = 
Moaul«[{xl. xZ« x3, x4, x5, x6, x7}, 
{xl, x2, x3} = TwoodJ.cElemoJa.tf 
x4 = t.Bncrt^[ ^3C'ii»J'^i'vePolynaiiii.alj -1; 

x5 = Hodule»ultplleatlon[KodulaHln?] {Pr-lmltlvePolynomi n1 ] [xl, xl] // Zaz'oPad[x4] / 
x6 = ModuloM>iltplication[«todMloIlins] [Prlnd.tlirePolyxioniJ.ol] [x2, x2] // ZezroPod [x4 ] ; 
BXod[2x6 -^xS, ModuloRlzi9] ] 

JUj.tomorxilulmaSlgma[MeduloRiJi{r_] [PrijaltlTrePolyzioiiilal^] :s 
Hodula[{xl, x2, x3, x4, xS, x6, x7, xB, x9}, 

xl - T*fOlVdlcExpozLBloxL[HodulaRiLn0] [Prl211i.tl1reP0lys.01ai.al] ; 
x2 = Hap [o[KoduloRJjiGr] [PrimltlirePolyiaonilal] , xl]; 

x3 = Map[ZeroPad[Z.en£rt:lk[PrlultlvePolyiiooiial] - 1] , Transpose [xl] I [3] ] ] ; 
{x2. 3e2} // Tranapos*] 

UIlltBRiag[ModuaoRlJa3_J [P3ri™itiiroPolyiiom±al_] : = 
Module [{xl, x2, x3, x4}, 

xl = AutonoxplilBiiiSlB>nB [ModuloRlng-] [PrlmltivaPolynomlal] ; 
x2 = xl // Trajoepoee // #[ [1]]6; 
x3 = x2 // 19od[2«< 4]G.// union: 
complaaent [x2, x3]] 

OnltBRlng[«] [{1, 0, O, 3, 2, 3}]; 

7-[Order_] [a_] [x_] r= (Neatl.lat [tr. x. Order- 2] // Apply [Plus, »]&) 

TraceRepr»sentatlon[Mod»loRlnff_] [Priaiiti-wePolynomial_] i = 
Module[{xl, x2, x3, x4, x5, x6, x7, xB . x9}, 

x2 = A\itaanorplilamSigi&a[Hod'uloKls9] IPrliiiltiirePolyTipnu.al] ; 
Map[(x3[Evaluata[#[[l]]]] i= Evaluate [#[ [2] ]]) t, x2] .- 
x4 = Hap[{«, 7-[xl] [x3] [«]}&, Tranepoee [x2 ] [[!]]] //Moa[#, ModuloRing] & 
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^->S Tupl«R«pr«««ntatlont»IoauloRijasr_] [Pirimiti^-BPoXynomiaJLI := ModMlertatr ac2 «3\ 

«JL = l.engrtb[ PrlwitivePolyaomialJi ^^e^ati, xz, x3), 

x2^ = ^NeetI.i8tjModuloMwltplicatlon[Mc>a^^ ^ P,.jjj^^^^p^3^y,^^^^^j ^^^^ ^^^^ ^^^^ 

x3 = Maipt zeros ad [xl -1], sc2]] 
This algorithm does not always yield thetcorrect value for the Minimum polynomial. 

W-Ln i mmnP oly[ModuloR^g_} lP^ijni.fcj.-.rePol3fnqmlal_] [rootPower- 1 fvar 1 t- 

Mod^.[{3cl. x2. x3, x4. x5. 3c6. x7 , xB, xlO. xll, xl2 , iLlB, xl^, xiS, tom>, 

xl = l.^3agrtH[ PximitlvePolynowial] -X; 

x2 = 2»1 _ 1; 
x3 = i»eBtt<ioc[Mod[# + #, x2]Bl, rootPowar, x2] r 
x4 = PoBitioji[x3, roofcPorwer] (»[[2, 1]] 
x6 = Talc«[x3, x4] ; 

3c9 = (xB, x7> // Tfansvosei 

xXO = Map[(Trai: + tam[#])ti, x6] / / ADI>ly [Timee. #] fc / / Expaad; 

niXe = com[n^] toin[ii^_] -> torn [n + mj r 

xix = xXO anxXB; 
xxa = TatoXeC Coef f Iclent [ xXi, t^ux, xX - i] , {i, o, xX}]y 
Map [ (torn [B»^XTia.tB[#[[X]] ]] : = EvaXuat© [# t [2] ] ] x91; 
ccan[xi_] t= EvaXuat«[ torn [Mod [n, x2] )]; Mod [ 
xX2 , KoduloRlng] ] 

Graef £eMetliod[aS'2inirq?oXy_] i ±= 
ModuXo[ 

{xl. x2odd, x2cnrexi, x3odd. x3e^exi, xAevea, x4edd, x5, xGeiren, x6odd< x7 xd 

= '^if ^1:: /= Floor)] // FXatten// TaXe[#. -xX]6; 

xze-w^a = jratoX«[{0. 1>, X, fxl + 1) /2 // PXoor}] // FXatten // Taleett. -xl]t; 

x3ev«a =. aF2mlnpoXyx2av«n // DropuaadingZeroBi 

x3odd = OF2iiiljQpoXy x2odd // DroplieadlngZez-oa; 
xAeves = PoXyiiaiiiiaXB«ultipXieatioii[4 ] [x3«v«ii, x3«v«ii] / 
x4odd a PoXyxioini.aXMuXt^d.plleatloa[4} [x3odd, x3odd] ; 
x5 = Bax[{t.engcli[x4e'vexa] , X.enErt:h [x4odd] } ] ; 
xGmnnaa. = ZeroPad [x5] Ix4©von] ; 

x6csdd = Zor-oPad[x5] Ix4odd3 
x7 = xCevezL — xEodd; 
xB = (x7 Sign [x7 [[!]]]) //Mod[«, 

Joiii[xB, {0>] // Partlt±e«a[», 2) &// Transpose [#][ [X] ] t ] 

Se<ittexicoGeafsirator[MoauXoRiiiff_] rPrimiti.vaPoXy_] [inltlalCoxidl rs©ar.eiiffth 1 
MoauXB[txX, x2State, xBDpdace, x4}, -j i a<-w_j 

xX = Mod [-Host [PrimitlvePoXy] , ModuXoRlzig^] ; 
x2ST^atie = XnltlaXCond ; 

x3C;pdate t= ModuXe[{}, x« = JCa^t [x2St:at:6] ; 

x2Stat::e = Joxzi [ {Hod [xSState . xX, HoduXoRi.ng] ) , l^rop rac2BT:ac«, -xn • x41 - 
{TabXd tstSnpdate, {±, X, 8eql.««istl») ] . x2Stata} } 

la±C±aXConditlcmB[Prijnlt:liraPoly_] := Module [{xX. x2. x3, x4, x5a, x6. x7^ 
XX = I,oiigrth.[PTrdLiiiibl-iroPoly] -X; " 

x2 = TBljXet{x3[i], D, X}, {i, X. Xl}]; x4 = Tatola [x3 [i] , (1, 1, xl}]; 
x5J = (ratolo ee Se«iuaiice[aoia[{x4}. x2]]) // Plattexi[#. xl - X] t; 
x6 = Mod[x5a[[l]] + 2x5J[[2]], «] ; 
x7 = Map[ ModtxSJt [2] ] + 2 #, 4]&, x5 JJ 
Joixi[x7, {x6}]] 

GoiaSe<iuejcice[Pir±ni±t:lvePoly_] [x&J.tco&dzrumbar_] := ISodule [ fxl, x2, x3 x4> 
xl = teiifftlilPrijiaXtivePoly] -1; " 
x2 = IjiitialConditioiiB [Primitl-vePoly] [ [XnltCoxxdNumbBr] ] ; 

x3 = 5ai3iienceGeiiQr'atior[4] [F7lmltltr»Poly] [x2] [2 -"^xX -1]; 

X4 = x3[[l]] /. (0 - 0, 1 - 0, 2 X, 3 -+ X) ] 

A-aeoeex'Z'elatlonSeixuQzicQ[GoldLSBi2uen.ce_] z = 
Module [{xX, x2Btiate< xSUpdate, x4}, 
xl GoXdSaizuance / . 0 -» -X ; - 
x2St:at:e = xl; 

x3X7pdat:e r= Module [{) , x4 = xX . x2StatB; xZState = RotateHigbfc [x2 State] i x41 j 
Table [xSupdate. (i. 1, Z.engrt:ii [GoXd3e<zuBZiCB] }] ] 
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8 



Moduie[{ica, yd, at3, } , 
xl = Xien0tli[PrlzadLtlv«Poly] - 1; 
3ca = znunJbKClBt. /. lsa.a^ -* 2''xl ->■ .1; 
Map[aol(S£l«q:u«xu;*[Px-liiiltlv»Poly] [«]&, x2] ] 

Module [{Id, x2, x3, x«}, 
acX = X.eziet:h [Fz-lmltilvePoly-] - T; 
sc2 = Tatolell. (1, 1^ 2-xl + 1>]; 
Kmp[oolds«Qu*nce[I>z-linltlirePoly] x2j ] 

Modula[{xl, jc2 State, x3 Update, x«}, 
xl s QoldBOQX/. 0 -» -X: 
x2Stats = (OoXdS*a2/. 0 ■* -X) i 

acSOpdaca i= Kod'aXB[(>, x4 = xX . xZState; x2State = RotateKlfrlit [x2Stat:e] ; x4] ; 
TabXe [xSTTOdate, (1, X, X.eiiertli [GoXdSe<zX] }] ] 

£xid[] 
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T: :-a.0a.gB = "TliiB le tbe symbol period" 



BT: iusage = 'Tbis i.B rlie usual produce 




B: :uaa0B = ^•riiXe dexiotea modluleibloD dLndLex Pi." 
<fr: s-txaage = '^[z.,e] la X.aurQzitia function" 

hFllt^sred: luaage s "Thla ia tihe fllCerBd gauaelaii pulse" 

PloaseAricrle: :uaaere = "This function takes eome time to calculate. Tbe Collowlng 

code will clraw a gprapb q£ tbe £uzic:ticnlPhaaeI'Unctlon[t_] 
phasexKSlncs =Tal>le [{t. PliaaeFunctlon[c T] } //N, {e, 0, I< ,1/40}]; 
ZilstPlot [pbasepolnte, Plot Joined -> Tru*] " 

Ptia.aeAxig'lePBet:: ruaa^B = "PhaaeAjaerleFaat [Z.] [t] speeda up tbe calculation 
of pJoaa^Anffle by ealculatlnff PbaseAn7la[ii] [t] wltli. I. numeric and t 
symbolic. Xt takes some time to calculate." 

E ; rusag-a = "S = Sla[»] " 

az lusa^e = = ■ ~ 

C::uaa0e="C = Cos [»] " 

M::uBa9e="H = 2*" " - 

Uodula.t:loziXxidex: z-amAga = "ModulatlonZndex = b- 

Z.BUJrexxtS: ludags = "LaUx-^ntSlX.] [t] = Sln[ <fr[l.,t + n *] ] / S- 
Z.aiurentC: luaage = 'l>aurentc [K] [t] 1b Z^aurentiB Cit,t " 
AlpbaKX; lUaa^e = -AlpbaKl [Xili] [X, 1] la taixarQnta OK.i. " 
LaurentUCt :uaasa = "I^iaa:emi:l-K[l-J [X] tjlvea tlie support o£ Cj^t " 



AKKTztltlalGtatsBetOp: r usage = 
''AI7KInltlalStataSatiT^[z.] [K] [TnltBltSeiz,AccumulatedPbaBe] sets vqc tibB aequence 
oC prior states oC Ak.b tbat tlie modulator went tbgougb fce ^et: to tlie 
conatelation point specified fay AccumulabedPhase wblcli. Is really Ao,o~ 

AKN::usaBQ= "AKH[I.] [K] [ {State. AccumulatedPbase} ] deflxisa a^.b 1a terms o£ Ag^s" 

UodulatlziffPulee: lusa^e = 
"The Pulse 1b asaumsd to Iiave tlie follovrln? structure Pulse t^-] [K] [t] " 

Rumberof Cur-tres: : usage = "The numbBo: o£ pulses used by tbe modulator" 



SawpllngrlPterval: iiiaage = 
- SanQsllnglnterval Is tbe Interval between sampled of tbe output o£ tbB modulator' 

TnltlalState: lUBB^e = "ZziltialState Is tb« eat of bits tl»at are assumed 
to be present before l>e {B-i#a_3, .-.)'' 

StartingQuadrant: :uBa?e = ^StartlngrQuadrant = Ao..iazid Is a nwnber' 



The start of Modulator Definitions 
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Moa.ulBt:or: zuaage = -Hodulac.o^rCXi] [Bi.t;S«Q.OpCB] assumes tihe t.ol.XavijD.'g 
} d«£aulc opclone stairclnEKiuaajraiit -» 0, InitlRlStat© -* Table [1,(1. 

X, 2 0} ] , S«KDli3i3lnt:.erval ^ T/32.sn3i]i2>erO£CU3rv'da -» A,HoduXat±3aaPuXB«-» 

Z.aurent:C . The Pulse B.&^\ixa^A to tuive the £ollowizi^ etnacture Pulde [I.] [K] [t] " 

Start of the Receiver Functions *^ ' 

riXtTulser ruBB^e n -The de£ault pulse and is caJLled fay PlltI»uloo [L] [K] [t] " 



Receiver: :usag« = " Rece±-«reir[Z.] [XxLEnit3eiz,C>ptB] asawneB tlie Collowln? 
de£ault options {SCaa:tli!isrQuadaraAt-».0, l±kitialStB.te-»-[l, 1, 1, 1, 1, 1, 1, 1, 
1, 1. 1. 1,1. 1, 1, 1«1. Xr 1. l},aaiiipliAgXxiterval -> 1/32 .iaaA\xla.t^±a.gVMXa^ -» 
FiltPulBe, SyncSanple -> O^NuinlMrOfCurves -> 2 - 

RoCol-»ei:Pxroi)oi:s luaasro = "aecoi-wer-Proiioi-[t.] I Sbax-tlsErQu.adz'azit, IxiltialBtati 
aaiBpl-ir>gTnter'«ral,ModuJ.&ti»gI>ttlBe. SyncSoxmile.ITutntoexOSCurveB, XnputseQ] 
Is called by Recalvei: a£ter all the optilorks ba-ve been zr^sol-red" 

[ * ' Pz-lvace • - ] 



2 ttBT 

Ideally we would and did deflne the effect of the convolution of h[t] by the formula below. However, this 
Mathemaiica gives on error. 

b«±i [t_] := Relaase[Hoaule[Ct) , J ^ «31t]j 

hK±lter-od[l»ulBoWielth^] [t_] := Modulo [(scl, x2, se3J, xl = 

. rr - h[tl-T:] , , I>iilaeWidth PulooWiath X 
(Mt»l. 40] ft) [Table [{tl. dr}. {tl. . . 1^}]]' 

x2 = li»t©xi>ol^tAoa[«3.] ; »2 [t] ] 



c := Cob [«] ; 
S := Six»t»l; 
a t= { // caiop) ; 
H := 2^-=^.- 

PhaaeAngle [!._] [t_] /rtS0t=O 
vhaseAzx^le [t_] / » t i I. T t = I 

PhaeeAnsrle [!._] [t_] i = 



x2 = Table [ 



e[{t2, 9 ^ Evaluate [xl] ditl), {t2, 0. L T. -5^^}]^ Interpolatioi* tai2] [t] ] 
We need to ensure that we always calculate PhascAnglc with t symbolic first to do the calculation only once 
Pb»BeaiiBleFaBt[I._] tt_] == Moohile [ {sel> , xlltt_] = PhaseAnglo [1.] [tt] ; jcl [t] ] ; 
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O'EI.-. t_] /jO<t<I.Xt= Pli«BeAag^leF»Bt [1-] [t] 

((rtl._, t_] /: 3ImT > b fc I.T := B - FliaaBAzigleFKBl: [I.] [t - I-T] 

We need to put ihis to avoid the function being extrapolated 

t_] /J I (0 < t < liT) tfi> J (2 LT > t i tT) := 0 
tawMts [e_] »= Sla[ ^^r[r., t + » T] ] / S 

AlphoKl tlX_] [K_, i_] /; (0 < i < IJ.) a*. (0 <= K < a"-*) i = 
KoaxiJLo[{icl, 3c2, 3c3, K»uin), xX := ^x2 = Mod[KIIuni, 2] / KNum = ^ ; 3t2j ; 

KNum.<> K.> x3 = Tabl*[3cl, {11, 0, - 1} ] / x3 fl]] ] 

■ AlpbaOCLri.] [X. li]) -i±. (l.±, X, X>]j Mlxi[xl]] 
I,aur»iitC[I._] rK_] tt_] /f 0 <= K< a*- : = 

l.aus:eiitS[I.] [0] [t] Y\ I^wewntS [X.] [11 + 1. AXpfcaKI [1.] [K. 11]] [t] 

i.i-=x 

The start of the modulator function 

AiqKTzLlt:laJLBfcat:eSe1::'Dp[Xi_] [K_l [XzilCBlcSBqt_, AccuinulaCed£>liaBe_] : = 
Module [{xX, x2, x3, xA, x5, acustoaae, Inl-tl^ltseq}, 
InlcblCB^q = inltBlbSea;' 

npaateseq 

tjoauleH), xX = acux>base - Svxn [InittitoBq [ [1 ] ] AipliaXI [L] [K, 1], {1, X, I. - 1) ] .- 
acupliase = acrupbiasa - Fir-et: [ln±i:bxt:seiz] ; IzilnblcsMQ □ Raab [ Inltl^ltsea] / xl] / 
Tal}Xe[ UpdateSeiz, (1, 1, I-aur-eatliK[li] [X] ) ] ] 



sum [atata [ [1 ♦ X] ] AlpliaKI [I.] [K, 1] , {1, X, I. - X) ] 



Modulator- [Z._] [BltSeii_, opca : = 

Module [ 

{xl, x2, ic3r tC, xB, x6. acaca, Aecuinalafced£>liaBe, aeii, AXMState, cuzTres, Pulse}. 
3cX = SaxnpllaffXxirerval /. {Opta} /- Qptlcmd [Kodulatior] ; 
atace = inltlalState /. {C^ta} /. opclosis [Modulator] i 
x3 = StartingQuBdrBiit / . {CSpCs} /. PpticwB [Modulator] ; 
x4 = Soaqplinglnteriral /. {Opta) /. Options [Modulator] ; 
X-vlXsq = MoaulatlnffSmlaa /. (Ol>ts> /. Opt xon^ [Modulator] ; 
Ciurv^aa =. (znunberOSCurveB V . {Opta) /. optloiis [Modulator] > - 1/ 
Acc\unuXatedPlia.Be = x3 ; 
ae(z = BltSaiz; 
Tablet 

AiagBtate[X] = jkMKInltialStateSetiJp[l,] [K] [state, AccumulatedPloaBe) , {K, 0, curves}]; 
x5 := Module [[}, state = Joiii.[ {Firat [seq] } . Drop[Btate, -1]]; 
AccumulatedX^aaa = AccuznulatQiiPlrLaee First [seq] / 
aeq = Raat [0e<z] ; 

Table [AKIVStata [K] = Jols[ [AKH [X.] [K] [ {state, AccumuXatedPbaae} ] } , 

Drop[AKNState[K] , -1]], {K, 0, Curves) ] 7 
x6lT_] = SmD[Sum[ (j)«»^"*-Wiri*il] inilae[Z.3 [K] [r + It], 

{1, 0, r.auroiitI.K[I.] [K] -1)], (K, 0, CurvOB}]; 
Table [xS [t] , (r, 0, T - xd, x4)]]; 
Tal>le[x5, {Mc, X, r.Biifftlii[BltBoq) ) ] / / riatto»] 

Options [Reeeiver] ts { Start IngQuadrant -» 0. 

InitlalState -» {X, 1, X, X, 1, 1, 1, 1. 1. 1, 1, 1, 1, 1. 1, X, 1, 1, X, 1), Sawplinglxitor' 
Modulatl&H£^Xse -► FiltPulao, SyncSaiqple -> O, Mumbero£Cur-vea -> 2); 
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Modtile C (acX, 3c2, ac3, x4. x5, xC}, 

xl = StartlnffQuaairant: /^.. {Outd} /. Options [Racsi vex-] ; 

x3 = Sanqplloglntez-ral /. {OpCs) /. ppcloxia [Racelver] ; 
xA = Kod«ilat:±ngPulBB / . {Opts} / . Options [RecedLver] ; 
x5 = ByzicSBiaplQ /. (Opte} /. Options [Recslwr] ; 
, x6 = Ruiob«rO£CurveB / . (Opts) / - OptlonB [Reealvar] • 
Kecelv<erFrep«z-[x.] [xl, x2, x3 . »cS, scG, lupvitSea] ] 

RscelverFrpper [!._} [ StartiiigrQuB.iraiit_, XndLtlalstata_, S«inplina'l»itBinrBi_, 

KoduladngPulsa . PyncBiUop3ke_, NumbezOfCixznreo^, 3:zwutS«9_] ; = 

Koctule[{xl, agn. Racei-vedseq. xxpectedvalue, ao(z# KeoelveKaxt, D, ^} . 
xl = T / SasvlljaglnteirvaX; 

RQcel-v-edSegc = Pa2rtltlon[ZnputSatz< xl] // Vxainapose // » [ [SyncSeunple -t-lJlK; 
ExpactedVaXuo = 

Modulatlngnilse [Xl] [0] [ (Z.aui:entUC[l.] [0] /2) T SyncSampla SastpllngXzitaxriral] ; 
X - Startln^QuAdz-ant; 
egn ~ Ot 

D = {}; 

seQ = Reeelvedsaa; 

Recel-veHexC : = ' ' 
Module [{xl, x3>, 

xl = { (-1)""" J* Fix-Bt[Beq]) / / Zsns 
. X£[Alia[jtl. - ExpectedValue 7 <= Mas [xl ♦ ExpactadValue ] , 

ID = iToijtt[ D, {l}]r J» = Si +X, D = aotn[ U, {-1)3; = ^ - 1 ] ; 

Beix = Rest [sect] J sgn = Med[0gB ■•■ 1, 2]]; 
Table [RecedLveHaxtr {i, 1, X.ezxstli.[Recel-vedSea] } ] ; 
ffi] 

• End[] 



EndPaelcaBn [ ] 



